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A model,  proposed  earlier,  is  modified  in  an  attempt  to  explain  a number  of  curious  behaviors 
of  corrosion-fatigue  crack  propagation  (CFCP).  The  behaviors  include  effects  of  load  ratio  R in  air 
and  salt  water  vs.  vacuum,  and  effects  of  loading  frequency  at  fixed  R in  these  environments. 
Assumptions  of  the  modelling  are  reviewed  in  detail  in  view  of  earlier  objections  to  them.  The 
ingredients  of  CFCP  per  this  model;  Poisson  contraction,  strain  hardening,  ligament  surface  attack/ 
annihilation,  and  stress  relaxation  are  developed  and  related  to  conditions  of  the  crack-tip  locale. 
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20.  Abstract  (Continued) 

In  the  modelling,  a parameter  G,  for  growth  rate  factor,  is  developed  solely  as  a function  of  the  form 
of  the  ordinary  or  of  the  cyclic  stress-strain  curve.  Previous  work  had  developed  a GJ  for  the  ordinary 
curve,  to  be  associated  with  the  surface  attack  effect  as  in  stress-corrosion  cracking,  and  one  G2  for  the 
cyclic  curve,  to  be  associated  with  the  stress  relaxation  effect  as  in  fatigue  crack  propagation  (FCP). 

A hybrid  G2m  is  developed,  combining  attributes  of  both,  which  seems  to  successfully  describe  the 
corrosion  induced  augmentation  of  G2vgParametric  curves  of  G 2/ 1 correspond  well  with  stage  II 
frequency -dependent  growth  in  CFCP.  Howevek,  alone  they  do  not  explain  the  frequency-wise  stage  II 
threshold  shift  nor  the  frequency-independent  air-environment  FCP  rate.  It  is  found  that  these  trends 
can  be  represented  by  loci  of  constant  plastic  strain  rate,  due  to  crack  loading  and  propagation,  relative 
to  the  surface  annihilation  rate.  Such  loci  are  determined  by  comparing  “G  maps”  with  strain  rate 
maps,  using  parametric  curves  of  equal  geometric-series  spacing.  Maps  of  this  sort  are  used  to  analyze 
about  a dozen  cases  of  CFCP  including  two  titanium  alloys,  one  aluminum  and  three  steels,  with  one 
of  the  steels  of  four  different  tempers.  Stress-strain  curves  of  the  low  strength  steels  are  processed  to 
remove  the  Liider  band  effect  to  facilitate  the  modelling.  The  scheme  for  data  organization  involves 
a representation  of  indexes  of  the  two  kinds  of  parametric  curves  fitting  the  data,  and  the  process  zone 
size  implied  by  the  fitting.  Comparison  of  estimated  process  zone  size  with  literature  data  of  micro- 
structural  and  fractographic  size  measurements  is  encouraging.  Model  predication  of  load  ratio  effects 
on  the  fatigue  crack  growth  threshold  is  in  good  correspondence  with  literature  data. 

For  the  reader  who  may  wish  to  defer  a detailed  study  of  this  report,  an  Appendix  provides  an 
illustrated  summary  of  major  points. 
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Curious  effects  in  CFCP:  of  load  ratio  R,  a)little 
effect  in  vacuum  but  b)  increased  growth  rate  and 
decreased  AK-threshold  in  air,  and  c)  increased 
stage  II  level  in  salt  water;  and  of  loading  fre- 
quency f,  d)  continuously  decreasing  stage  II 
onset;  e)  stepwise  decrease,  and  f)  no  change. 

Various  plastic  flow  terms  as  determined  from  the 
true  stress-strain  envelopes  obtained  from  the 
flow  property  test. 

A straight  line  fit  of  a logarithmic  plot  of  true 
stress  vs  strain,  of  slope  n,  provides  an  equa- 
tion for  use  in  the  LUder  band  region  of  the 
medium  strength  steels  of  this  study. 

The  "de-LUderi zation"  procedure  includes  a)  stress- 
strain  curve  for  X-65  steel  with  upper  yield  point 
in  first  (quarter)  cycle;  which  b)  converted  dir- 
ectly to  growth  rate  factor  maps  shows  singular- 
ity; while  c)  converted  via  equation  from  Fig.  3 
is  continuous;  whence  d)  combined  front  from 
b)  and  rear  from  c)  provides  G(N)  map  used. 

Combined  growth  and  strain  rate  maps  for  X-65 
steel:  a)  Growth  Rate  Factor,  N-map,  from  Fig.  1+ ; 
b)  parametric  set  of  lines  of  constant  plastic 
strain  rate,  M map,  in  coordinates  of  a);  super- 
posed in  c)  whence  segmented  curves  of  constant 
strain  rate/surface  attack  rate  N-M  map  are  drawn. 

The  minor  effect  of  load  ratio  R on  FCP  rate  in 
vacuum  is  compared  to  the  m-induced  G^  growth 
prediction  of  the  present  scheme  for  the  Ti-6A1- 
1+V  noted. 

When  exercised  in  air,  FCP  rates  of  Ti-6A1-Uv 
are  markedly  increased  and  the  AK-threshold 
decreased  with  R;  the  growth  pattern  is  matched 
by  curves  of  constant  ratio  of  total  plastic 
strain  As  to  surface  attack  Arp  in  each  load 
cycle,  N-M  = 9 set.  At  the  highest  R,  these 
curves  are  poorly  defined  by  a line  of  constant 
growth  rate  factor,  N =-1.5  matches  here.  Note 
that  pertinent  data  on  this  and  other  case  studies 
is  given  in  Tables  indentified  with  Figure  num- 
bers. 
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Salt-water  CFCP  data  of  Pettit  et  al  for  a Ti-6A1- 
UV  is  bounded  by  frequency  independent  AEp/Arrp 
lines  of  N-M  indices  shown  for  Stage  II  onset  as 
an  upper  limit  of  growth  rate  (N-M  - 15)  and  air 
environment  as  a lower  limit  (N-M  = 9 as  in 
Fig.  T results)  whereas  crossover  is  at  constant 
Arip  for  the  cycle  (N  set). 

Salt  water  CFCP  data  of  Meyn  for  a Ti-8Al-lMo-lV 
is  bounded  in  the  manner  of  Fig.  8.  The  lower 
bound  G2~only  curve  of  the  predictive  map  is  of 
the  same  trend  as  Meyn ' s FCP  rate  in  vacuum  but 
fails  to  anticipate  the  trend  at  low  growth 
AK  levels. 


Dry  argon  and  salt  water  FCP  data  of  Bucci  for  the 
Ti-8Al-lMo-lV  of  Meyn  shows  the  argon  data  fitted 
by  about  the  same  Aep/Arp  (N-M  - 10.5)  as  the  air- 
environment  data  of  Meyn.  The  salt  water  CFCP 
trends  are  fitted  by  constant  Arp  per  cycle, 

N - 10,  for  all  thru  R levels  with  no  apparent 
Aep/Arp  upper  limit. 

Dry  argon  FCP  data  of  Unangst  and  Wei  is  poorly 
matched  by  G2  curves,  a failure  unexpected  in 
view  of  accurate  model  predictions  of  FCP  for 
many  other  aluminum  alloys. 

Fresh  (distilled)  water  and  air  data  of  NRL  on 
4 3 U 0 steel  of  four  tempers  is  shown  in  the  next 
four  figures.  Here  the  hardest  condition,  205°C 
tempering,  appears  nicely  bounded  N-M  = 19  for 
upper  frequency-independent  limit,  13.5  for  the 
lower  one,  with  various  constant-N  lines  fitting 
stage  II  data  in  transition  there 

G-map  vs  CFCP  data  for  U3^0  steel  as  in  Fig.  12 
but  for  tempering  at  315°C.  FCP  data  of  Gallagher 
is  much  lower  than  NRL  data  for  air  but  not  quite 
as  low  as  present  prediction  for  m-induced  growth, 
the  G2~only  line. 

G-map  vs  CFCP  data  for  1+3^0  steel  as  in  Fig.  12 
but  for  tempering  at  h2'J°C.  A marked  decrease 
in  the  level  of  stage  II  growth  rate  is  evident 
here  . 

G-map  vs  CFCP  data  for  ^3^0  steel  as  in  Fig.  12 
but  for  tempering  at  538°C.  The  G-map  has  been 
prepared  by  the  "de-Ltlder  i zing"  process  to 
remove  the  upper  yield  point  effect. 
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. 1 6 : Two  data  sources  on  a 5-Ni  steel  are  compared  with 

6l  a G-map  for  a low  R-value  - 0.1.  In  both  cases 

the  lower  bound  is  matched  by  curves  of  G ^ only, 
indicating  a negligible  effect  of  the  air  environ- 
ment. The  upper  CFCP  rate  bound  and  transitional 
curves  are  as  in  preceding  figures. 

. 17:  CFCP  data  on  X-65  steel  is  roughly  matched  by 

6 3 G-maps  developed  in  Figs.  3,^,  & 5-  However, 

this  strength  level  is  thought  to  violate  the 
lower  limit  of  applicability  of  this  model  as 
presently  constituted  without  reckoning  for 
crack  closure  effects,  and/or  microstructural 
vs  dT  size  effects. 

. 18  : Data  organization  is  proposed  in  terms  of  this  log 
6 5 plot  of  Ar<p  vs  cycle  duration  x,  plus  the  ASp/Ar.p 

values  (N-M  indexes)  defining  frequency  indepen- 
dent limits  of  Table  II.  Here  the  one  second 
intercept  is  compared  with  Vg  estimates  in  stress 
corrosion  cracking  experiments.  A 1:1  slope 
denotes  constant  Vg. 

. 19:  The  model,  as  a result  of  T-C  envelope  centroid 

67  neutralization,  predicts  the  trends  shown  for 

the  AKrpjj  limit  vs  R,  which  seems  fairly  consistent 
with  extensive  NEL-data  of  Pook. 

. 20:  Comparison  of  calculated  dT  process  zone  size 

69  with  measured  fr act ographi c /mi cro struct ur al 
features,  from  sources  noted  (see  Table  III) 
shows  a reasonable  correspondence. 


ORGANIZATIONAL  SCHEME  FOR  CORROSION-FATIGUE 
CRACK  PROPAGATION  DATA 


I.  INTRODUCTION 

This  paper  continues  authors'  pursuit  of  physical  models 
by  which  various  manifestations  of  subcritical  crack  propa- 
gation can  be  predicted  from  ordinary  and  cyclic  stress- 
strain  properties  of  the  cracked  material.  The  materials 
for  which  this  now  appears  feasible  are  generally  high  to 
medium-strength  wrought  alloys,  whose  fracture  instability 
involves  ductile  tearing,  albeit  localized,  as  distinct  from 
cleavage.  A central  assumption  of  this  modelling  is  the 
necessity  and  likely  sufficiency  of  condition  for  tensile 
instability  in  the  crack  locale  for  fracture  instability. 

The  original  proposal  (196U)  [l]  stemmed  from  experimental 

correlations  of  the  strain  for  simple  tensile  instability, 
as  influenced  by  temperature  and  strain  rate,  with  the 
plane-strain  fracture  toughness,  Kjc , as  similarly  influ- 
enced. Some  years  later  [2]  it  was  noticed  that  crack  vel- 
ocity under  conditions  for  stress-corrosion  cracking  (SCC) 
could  be  related  to  the  shape  of  the  entire  stress-strain 
curve.  With  strain  near  the  crack  tip  assumed  proportional 
to  the  stress  intensity  factor,  this  corresponds  to  an 
applied  K less  than  Kjc , i.e.,  the  regime  of  subcritical 
crack  propagation.  Here  the  source  of  tensile  instability 
was  viewed  as  a surface  annihilation  - or  corrosion-rate 
acting  upon  critical  crack-tip  ligaments.  Since  corrosion 
is  a stable  process,  so  correspondingly  is  the  crack  growth 
induced  to  compensate  for  the  eroding  sectional  area  of  the 
ligament.  In  a sense,  then,  it  became  unnecessary  to 
specify  tensile  instability  as  a failure  condition.  Crack 
growth  could  be  perceived  as  a mechanism,  other  than  rising 
external  load,  by  which  the  load  upon  crack  tip  elements  is 
mai n t a i ne d constant . Later,  at  the  suggestion  of  Landes  and 
Wei  [3],  the  stress  relaxation  during  the  dwell  time  of  a 
fatigue  cycle  was  recognized  as  another  source  of  load 
diminution.  When  related  to  cyclic  stress-strain  proper- 
ties, this  provided  a one-parameter  fit  of  a number  of  sets 
of  fatigue  crack  propagation  (FCP)  data  [ U ] . This  was 
followed  by  two  NRL  studies  [5,6],  attempting  to  combine 
influences  of  ligament  corrosion,  as  a surface  effect,  with 
its  stress  relaxation,  as  a volume  effect,  to  explain 
corrosion  fatigue  crack  propagation  (CFCP). 


Note:  Manuscript  submitted  June  13,  1977. 
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Only  partial  success  was  enjoyed  from  these  endeavors. 
One  difficulty,  carried  over  from  the  stress-corrosion- 
cracking model,  was  a failure  to  predict  the  correct  trend 
in  the  AK-threshold . Whereas  the  model  related  the 
threshold  directly  to  the  yield  strength  level,  the  data 
indicates  a strong  inverse  trend,  except  at  very  high 
strength  levels.  Another  deficiency  was  pparent  in  CFCP. 
Here  the  K-threshold  of  rapid,  environmentally-induced 
Stage  II  growth  tends  to  increase  with  increased  cycle 
duration.  The  model  predicted  no  change. 

These  shortcomings,  as  well  as  concern  with  the  sim- 
plistic analysis  and  assumptions,  led  to  its  less  than 
enthusiastic  acceptance  by  the  technical  community.  In  an 
attempt  to  remedy  these  shortcomings,  this  paper  reports 
results  of  a careful  examination  of  some  dozen  sets  of 
corrosion  fatigue  crack  propagation  data.  These  exemplify 
each  of  the  currently  known  anomolous  behaviors  of  CFCP. 

In  each  case,  we  have  performed  mechanical  stress-strain 
tests  of  the  actual,  or  essentially  identical,  materials. 
Between  these  end  points,  fracture  vs  flow  data  sets,  we 
have  permuted  what  combinations  of  model  parameters  it 
seemed  permissible  to  vary.  The  permutations  involved 
calculation  and  graphing,  made  possible  by  digital  computer 
and  graphics  programming.  It  is  thought  that  some  improve- 
ments in  modelling  have  been  developed  which  are  worth 
reporting  at  this  time. 

2.  ANOMALOUS  BEHAVIORS  IN  CFCP 

It  may  be  helpful  at  the  outset  to  list  each  character- 
istic of  CFCP  which  has  been  addressed,  citing  here  the 
data  source  and  material  taken  as  illustrative.  Such 
initial  summary  of  behaviors  follows  the  practice  of 
Speidel  [7]  and  of  Dawson  and  Pelloux  [8].  Figure  1 depicts 
those  which  are  treated  in  this  paper.  We  follow  the  Paris 
convention  [9]  of  displaying  FCP  data  as  a full  logarithmic 

plot  of  the  FCP  rate  = aN , vs  the  opening  mode  stress 

AN  " 

intensity  factor  range  AKj , simply  called  AK. 

Consider  first  effects  of  the  load  ratio  R on  CFCP  at 
constant  loading  frequency.  The  base  of  reference.  Figure 
la,  is  fatigue  crack  propagation  in  vacuum,  for  which  we 
employ  the  results  of  Irving  and  Beevers  [10]  on  a marten- 
sitic Ti-6Al-UV  alloy.  Generally  the  FCP  rate  is  lower  in 
vacuum  than  in  air,  particularly  at  low  AK  levels.  Despite 
the  indication  of  an  environmental  effect  of  air,  the  FCP 
rate  in  air  is  insensitive  to  loading  frequency,  atypical 
of  that  in  liquid  corrode  ts . Irving  and  Beevers  find  the 
threshold  for  onset  of  FCP,  AK-j^,  relatively  insensitive  to 
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load  ratio  R above  about  0.3.  In  one  ease  of  a low  alloy 
Q&T  steel  Cooke,  Irving,  Booth  and  Beevers  [ll]  observe 
practically  nil  sensitivity  for  all  R's.  In  normal  air 
environment,  however,  AKrp^(R)  drops  off  with  (increasing) 

R to  a limit  of  about  one-third  that  for  R = 0.  The  example 
of  this.  Figure  lb,  is  a Ti-6A1-^V  used  by  Fitzgerald  and 
Wei  [12].  In  mid-range  AK  levels,  the  FCP  rates  tend  to 
merge  again,  prior  to  AKmax  reaching  the  Kc  instability 
limit.  This  is  observed  for  the  Ti-6A1-W,  as  well  as  for 
A533  B steel  (not  used  here)  of  Paris  et  al  [ 1 3 ] - However, 
in  Ti-8Al-lMo-lV , Fig.  lc , Bucci  [lU]  observes  a strong 
leftward  translation  of  argon-environment  growth  curves  with 
increased  R.  Correspondingly,  his  salt  water  growth  curve 
shifts  even  more  to  the  left  and  upwards  as  well. 

Regarding  effects  of  frequency  at  fixed  R,  Pettit, 

Krupp,  Ryder  and  Hoeppner  [ 1 5 ] » Fig.  le , show  a typically 
increased  CFCP  rate  in  Stage  II  growth  with  cycle  duration 
on  a Ti-6Al-kV  in  salt  water.  However,  the  AK  for  Stage  II 
growth  acceleration  increases  with  cycle  duration,  leading 
to  curious  crossovers  in  CFCP  rates  in  the  transition  to 
the  normal  (air)  growth  rate  curve.  This  effect  is  clearly 
documented  by  Dawson  and  Pelloux  [8]  for  Ti-6Al-6V-2Sn . 

In  a high  strength  hsho  steel,  little  sweep-back  in  Stage  II 
threshold  is  observed.  Fig.  Id,  whereas  in  a line-pipe 
steel,  even  more  is  observed  by  Vosikovsky  [Uo],  Fig.  If. 

Dawson  and  Pelloux  showed  CFCP  for  square  wave  loading 
in  salt  water  about  3X  faster  than  for  sine  wave.  However 
Vosikovsky  has  one  example  in  X65  line-pipe  steel  showing 
a CFCP  retardation  by  1/3  for  square  vs  triangular  loading 
patterns.  Surely  there  are  other  anomolous  behaviors  of 
CFCP.  But  even  with  those  noted,  small  wonder  that  con- 
fusion is  rife  amongst  engineers  seeking  to  use  the  avail- 
able data  for  safe-life  prediction. 

3.  BRIEF  OF  MODELLING  ASSUMPTIONS 

In  view  of  aforenoted  criticism  of  model  assumptions, 
a word  here  on  rationale.  This,  like  any  physical  model 
for  fracture  in  an  active  environment,  requires  a failure 
criterion  in  the  context  of  a crack  tip  stress  and/or  strain 
field  analysis.  To  fit  any  such  model  to  CFCP  data  will 
involve  a size  parameter,  implying  something  of  the  alloy 
microstructure.  Further,  account  for  corrosion  effects  will 
involve  some  sort  of  surface  attack,  implying  aspects  of 
corrosion  chemistry.  To  combine  all  of  these  elements  into 
single,  workable  fracture  model  requires  major  simplifi- 
cations/approximations of  each.  Such  simplifications  often 
prove  unacceptable  to  experts  in  the  field  involved,  leading 
to  rejection  of  the  entire  scheme.  We  are  not  astute 
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enough  to  avoid  this  dilemma,  but  ask  that  suspected  assump- 
tions be  granted  trial  as  working  hypotheses,  whose 
plausability  would  be  enhanced  by  success  of  data  correla- 
tions which  employ  them. 

A . Failure  Criterion/Plasticity  Model 

Regarding  the  failure  criterion,  we  hypothesize,  as 
Orowan  years  ago  [ 1 6 ] , that  local  tensile  instability, 
generally  conceded  to  be  a necessary  condition  in  the  path 
to  rupture  of  a tensile  specimen,  is  also  a very-nearly 
sufficient  one  for  the  crack  tip  locale.  How  nearly  depends 
of  course,  on  the  effective  compliance  of  links  connecting 
the  critical  element  to  the  interior  of  the  cracked  body. 
The  contention  that  reaching  the  maximum  load  point  is 
sufficient  was  based  on  the  experimental  evidence  of  direct 
correlation  of  the  plane  strain  fracture  toughness  Kjc  with 
the  plastic-flow  strain  for  maximum  load,  as  estimated  from 
the  strain  hardening  exponent  n.  Both  Kjc  and  n were 
observed  to  respond  similarly  to  variations  in  temperature 
and  straining  rate.  A size  parameter  is  derived  from  a cal- 
culation of  the  distance  from  the  mathematical  crack  tip 
origin  at  which  the  tensile  y-strain  level  recedes  to  the 
value  of  n.  With  Kpc  and  n directly  proportional,  a 
y-strain  directly  proportional  to  K (and  independent  of 
yield  strength  and  strain  hardening  characteristics  of  the 
material)  was  seen  as  indicating  a constant  size  parameter. 
Since  microstructural  size  parameters  should  be  unaffected 
by  reduced  temperature  and  high  strain  rate,  the  test- 
variables  employed,  this  (constancy)  seemed  a desirable 
attribute.  Hence,  the  simple  elastic  field  analogue  for 
crack  tip  strain  was  adopted 


ey  = KI/v/27rr  E 


(1) 


where  Sy  is  strain  normal  to  the  crack  plane  at  distance  r 
forward  of  the  crack  tip,  and  E is  Young's  modulus.  Since 

”jc,  the  size  parameter  was  calculated  as 


Ey  = n when  Kj 


<lx  = Kj2c/27rE2n2 


(2) 


B.  Microstructural  Size  Parameter 


Significance  of  the  size  parameter  dip , called  a 
process  zone  size,  could  not  be  assessed  from  evidence  then 
at  hand.  The  strain  hardening  exponent  n is  not  always 
equal  to  the  strain  for  tensile  instability.  Crack  tip 
triaxiality,  as  discussed  by  Williams  and  Turner  [ 1 T ] , can 
require  a greater  strain.  More  too  would  be  needed  if 
failure  required  attainment  of  a strain  intermediate  that 
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for  maximum  load  and  that  for  rupture.  Larger  critical 
strains,  at  given  Kpc,  are  found  only  closer  in  toward  the 
crack  tip,  hence  smaller  drp  . Conversely,  a stronger  than 
inverse-half-power  strain  singularity,  indicated  by  several 
analyses  [l8,  19,  20 ] , would  provide,  for  given  Kpc>  strains 
reaching  the  increased  "failure"  strain  at  greater  distances 
from  the  crack  tip.  Hence  a judicious  combination  of  these 
assumptions  could  provide  almost  any  process-zone  size 
needed . 

As  it  turns  out,  the  size  needed  well  may  be  simply 
that  indicated  by  Equation  (2),  as  a number  of  investi- 
gations, to  be  discussed  in  Section  22,  have  shown  fracto- 
graphic  dimple  size  in  various  alloys  and/or  largest- 
inclusion  spacing  in  steels  to  correlate  reasonably  well 
with  the  dip  Df  Equation  2.  Further,  careful  experimental 
work  of  Liu  and  his  students  [21]  has  invariably  indicated 
the  inverse-half-power  strain  singularity.  We  conclude 
from  this  that  the  simple  e las t i c -anal ogue  strain  singul- 
arity of  Eq.  (l)  is  the  better  choice  for  these  problems, 
which  is  thus  employed  here. 

C . Environmental  Attack 

Finally  in  the  list  of  model  ingredients,  the 
environmental  influence  parameter  of  surface  attack.  We 
would  describe  this  broadly  as  any  environmental  activity 
about  the  crack  tip  dip-element  which  has  the  effect  of 
gradually  annihilating  the  load-bearing  capability  of 
material  at  its  surface.  The  most  obvious  means  of  removing 
material  was  thought  to  be  dissolution,  which  was  then 
proposed  in  Ref.  [2].  Corrosion  experts  object  to  this 
degree  of  specificity,  pointing  out  that  other  processes 
are  known  to  influence  SCC  and  CFCP.  Thus  making  the 
cracked  specimen  either  anodic  or  cathodic  with  respect  to 
its  neutral  condition  is  known  to  accelerate  the  SCC,  hence 
the  CFCP  rate.  Added  stress  from  pressurized  voids  [22], 
eased  dislocation  glide  by  matrix  solute  interaction  [23], 
micro-hydride  swelling,  flaking  [2U]  are  all  possible 
mechanisms.  Unfortunately  this  model  cannot  discriminate 
among  annihilation  mechanisms.  However,  it  does  allow 
quantification  of  the  effect  of  such  mechanisms  in  terms 
of  the  effective  total  depth  of  surface  annihilation  and, 
if  normalized  with  respect  to  cycle  duration,  of  its  rate 
of  attack . 

I*.  SOURCES  OF  LIMITED  TENSILE  INSTABILITY 

Assuming  that  the  condition  of  constant-load  mainte- 
nance could  closely  approach  a sufficient  condition  for 
crack  growth,  the  following,  previously  published. 
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derivation  is  repeated  here  [25].  The  analysis  pertains  to 
any  unrestrained  tensile  specimen,  such  as  those  of  size  irp 
thought  to  regulate  the  rate  of  crack  propagation.  Tiie 
question  is:  what  amount  of  straining;  hence  strain  hard- 
ening; liencq,  in  the  context  of  the  crack  tip  strain  field 
at  constant  K , crack  extension,  is  required  to  maintain  a 
constant  load,  hence  stability,  in  a given  ligament?  Such 
load  P supported  over  area  A defines  a stress 

o = P/A  ( 3 ) 

o r 

P = 75A 

Differentiating  then 

dP  = u dA  + A do  ( k ) 

At  constant  load  dP  = 0 whence 
- o dA  = A do 

The  total  differentials  of  area  A and  of  stress  5 can 
be  separated  into  constituent  partials.  The  familiar  ones 
are  Poisson  contraction  and  strain  hardening  respectively, 
ingredients  of  the  ordinary  mechanical  tensile  instability 
phenomenon.  Less  familiar  partials  are  the  dT-ligament 
surface  corrosion/annihilation,  and  stress  relaxation. 
Consider  now  each  influence  in  detail. 

A . Poisson  Contraction 

Lateral  contraction,  hence  areal  diminution,  due 
to  elongation  of  a ligament  is  characterized  by  the  Poisson 
ratio 

v = ~ed/eg  ( 6 ) 

where  is  diametral  strain,  e £ is  longitudinal  strain  of 
the  ligament.  For  simplicity,  the  subscript  £ , designating 
a longitudinal  property,  is  deleted  in  further  development. 

For  a ligament  of  radius  rT , (or  diameter  dip)  A = irr,  ‘ 

or  , 

dA  = 2ffrT  drT  = 27rr.p  ded 

Hence,  the  area  partial  due  to  Poisson  contraction  is 
defined  as 

dA,,  = 2mrp  dcd  = ~vv  dj  de  ^2  ( ' 
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which  for  plastic  flow  £p,  where  v = 0.5  becomes 

3A„  = -fd$dep  (9) 

B .  Strain  Hardening 

It  appears  from  these  modelling  attempts  that  only 
the  plastic  part  of  the  measured  total  (elastic  plus 
plastic)  strain  hardening  rate  is  effective  as  stabilizing 
the  three  "inducers"  of  crack  growth:  the  surface  attack 
rate,  the  stress  relaxation,  and  the  Poisson  contraction. 

In  the  case  of  stress  relaxation,  the  reason  for  this  is 
clearly  the  plastic  (only)  strain  rate  sensitivity  of  the 
flow  stress,  which  is  equated  to  the  stress  relaxation. 

For  the  other  inducers,  apparently  only  the  permanent,  non- 
reversible  portion  of  strain  hardening  has  a stabilization 
attribute.  Hence  for  the  plastic  strain  hardening,  the 
stress  partial  is 


O</0  = dfp  = 0V  dtp 

(It  p 

C . Ligament  Surface  Annihilation 

Effects  of  corrosion  rate  on  tensile  instability  are 
usually  negligible  for  ordinary  macro-specimen  sizes. 

However  for  micro-size  ligaments  at  the  crack  tip,  the 
effect  can  be  substantial.  From  Equation  (7)  above  let 
drip  = -Vgdt  , where  Vg  is  the  annihilation  rate  at  the  dp 
ligament  surface,  whence  this  partial  of  A is 

DAS  = -jrdT  Vs  dt  (11) 

The  surface  annihilation  rate  Vg  need  not  be  con- 
stant, but,  as  in  ordinary  corrosion,  typically  diminishes 
with  time,  as  indicated  later.  It  is  convenient  to  des- 
ignate a symbol  for  the  total  surface  intrusion  during  the 
entire  cycle  period:  Arp  = Vg At  . 

D . Stress  Relaxation 

The  remaining  partial  of  5 is  that  due  to  stress 
relaxation.  The  effect  in  structural  alloys  at  ordinary 
temperatures  may  be  closely  associated  with  the  plastic- 
flow  stress  sensitivity  to  rate  of  plastic  ( only ) de f ormat i on. 
For  present  purposes  it  is  important  to  have  a realistic, 
yet  hopefully  simple  characterization  of  this  effect,  parti- 
cularly of  how  it  is  related  to  the  position  in  the  stress- 
strain  cycle  in  which  it  occurs.  Since  there  is  a consid- 
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erable  literature  on  the  strain  rate  sensitivity  which 
could  be  brought  to  bear  on  this  question,  the  basis  of 
conversion  should  be  noted.  Derivations  in  the  literature 
provided  by  Guiu  and  Pratt  [26],  and  later  by  Hart  [27], 
are  abstracted  here  for  convenience  of  the  reader,  using 
symbols  of  Guiu  and  Pratt. 

In  a tensile  test,  the  crosshead  speed  y is  matched 
by  the  elastic  plug  plastic  lengthwise  deformation  rate  of 
the  specimen  4 + £.  plus  the  elastic  deformation  of  the 

machine  z,  particulErly  of  highly  stressed  elements  of  the 
load  train  such  as  the  load  cell  and  grips 

y = C + t + z (12 

p e 

With  crosshead  arrest,  y = 0 whence 

Kp=-(Ve+i>  (13: 

The  machine  deformation  is  inversely  proportional  to  its 
sti f f nes  s 

s-£  du: 

dz 

assumed  constant  for  small  changes  in  load  P.  Converting 
length  changes  l to  (engineering)  strain  rates  e,  Eq.  (13) 
becomes 

,15: 

where  lQ  and  a0  are  initial  length  and  section  area  of  the 
specimen.  The  elastic  deformation  rate  can  be  expressed 

ee  = d/E  ( 16  ' 

where  E is  Young's  modulus.  Combining  the  specimen  plus 
machine  compliance  into  total  compliance  C^,* 


* Designation  here  of  compliance  by  C,  as  used  by  Hart, 
rather  than  by  M of  Guiu  and  Pratt,  our  usual  practice. 


whence 


p 


-CT<i 


(18) 


Equation  (l8)  says  that  a plot  of  o vs  o,  is  equivalent  to 
a plot  of  a vs  e . Hence  the  stress  relaxation  test  pro- 
vides direct  measure  of  o(e). 


Whether  determined  by  stress  relaxation  vs  time  or 
by  stress  vs  strain  rate  measurements,  the  character  of 
this  behavior  is  important  here.  Guiu  and  Pratt  [26]  assert 
that  the  absolute  sensitivity  of  flow  stress  to  relative 
(log)  strain  rate,  — — • is  constant  over  widely  ranging 
ti me/strain  rates.  J?®6.  Campbell  [28]  normally  displays 
his  strain  rate  and  strain  rate  history  effects  in  this 
manner.  This  is  consistent  with  the  observation  [29]  that 
the  absolute  sensitivity  in  ferrous  alloys  tends  to  be  con- 
stant for  varied  strength  levels.  However,  when  one  looks 
at  a given  alloy  over  a range  of  flow  stress 
duced  by  straining,  the  relative  sens'tivity 
appears  more  nearly  constant.  Thus  plots  of 
true  flow  stress  vs  plastic  strain  rate  for  parametric 
variation  of  strain  level  tend  to  be  more  nearly  straight 
and  parallel  to  each  other  in  full  logarithmic  display, 
than  semilogarithmic . Hart's  work  supports  the  latter  con- 
tention. It  is  commonly  used  in  superplasticity  calcula- 
tions where  very  large  deformations  must  be  estimated.  Our 
own  observations  support  it.  Thus,  we  have  elected  the 
latter  option  in  this  model  development,  using  the  super- 
plasticity literature  convention  of  designating  the  exponent 
of  sensitivity  as  the  constant  m, 

_ (1  l»g(o)  ( 19  ) 

d log  (ep) 

or 

m = d?/a.,  ■■ 

dep/ep  (20) 

If  indeed  the  strain  rate  sensitivity  is  of  the  kind 
characterized  by  Equation  (20),  it  is  possible  to  measure 
it  from  a log  plot  of  stress  vs  the  relaxation  time.  Assume 
here,  as  implicit  in  Guiu  and  Pratt,  and  Hart,  a path-  or 
h i s tory- i ndependent  mechanical  equation  of  state.  Consider 
then  a path  in  stress/plastic-strain-rate  "snace"  to  a 
fixed  plastic  strain  at  two  different  but  constant  strain 
rates  + dc  in  correspondingly,  two  different 

times  tj  >t<j*  The  stress  differential  between  these  end 
points  is  given  by  rearrangement  of  Eq.  ( 20 ) . 


dlogt 
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(21) 


_ llf  ^2  ^1 

do  ~ mo  — = mo  : 

r f, 

However,  since  the  terminal  (plastic)  strain  e is  the  same 


e 


whence 


do 

a 


m 


c/t->  e/tj 
6/t| 


do  = 


mo  dt 
^2 


-in 


_dt 

t2 


(22) 


(23) 

(2M 


til  (o)  = 111  ill  (t) 


(25) 


Thus  the  slope  of  the  log  plot,  noted  above,  provides  a 
direct  measure  of  (-)m.  This  is  much  more  convenient  and 
experimentally  more  accurate  than  measuring  and  plotting  the 
slope  of  the  d vs  time  data  plot.  We  find  good  agreement 
between  the  two  methods  of  analysing  stress  relaxation  data, 
and  where  available,  with  direct  strain  rate  sensitivity 
measurement  s . 

It  should  be  noted  that  the  above  derivations  assume 
that  relaxation  of  plartic  strain  as  well  as  the  strain 
hardening  associated  with  it  are  negligible;  however,  cor- 
rections for  both  effects  are  simple  enough.  The  strain 
hardening  correction  is  handled  graphically  in  our  experi- 
mental procedure  by  projecting  back  to  the  elastic  line 
with  an  appropriate  slope  reflecting  the  strain  hardening 
rate. 


This  result  permits  the  fourth  stress  partial  to 
be  designated  as  s imply 


On. 


_ dt 
mo  — 


(26) 


where 


t^/t^  is  the  ratio  of  the  hold  or  dwell  time 


(27) 

(28) 

tj|  to  the 
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loading  time  of  the  loading  waveform. 

E . All  Sources  Together 

Combining  all  partials  into  the  "load-maintenance" 
Equation  (5)  gives 

-d(0A,,  + 0AS)  = A (day  + 0om)  ( 29  ) 


Inserting  the  above  defined  expressions  rives 


(f  4 d^P  + *dTVs  dl)  = A (°P  deP  ~ m(J  y) 

O 

We  substitute  iri  for  the  s pec imen / 1 i gamen t area  A,  assumed 
circular.  It  is  also  useful  to  rearrange  to  separate  the 
"stabilizing"  strain  differential  required  to  maintain  the 
constant  load  condition,  since  this  will  be  traded  for 
crack  growth  with  the  crack  tip  plasticity  model.  With  this 


dtp  = 


ivs 

dT 


dt  + 


- ^ - > 


-l 


(31) 


This  equation,  while  developed  for  the  case  of  a tiny  liga- 
ment of  radius  r^,  is  applicable  to  any  tensile  specimen, 
wherever  subjected  to  tension  strain,  corrosion,  strain 
hardening,  and  stress  relaxation.  It  assumes  no  attributes 
of  the  crack  tip  locale.  But  this  is  now  required. 


5.  CRACK-GROWTH  VS  STABILIZING-STRAIN  TRADEOFF 


It  is  necessary  now  to  ask  how  strain  is  applied  to 
microstructural  tensile  elements  encroached  by  an  advancing 
crack  which  is  periodically  unloaded,  then  reloaded.  For 
reasons  noted  earlier,  we  employ  the  simple  analogue  of  the 
linear  elastic  crack-tip  stress  field.  Repeating  Equation 
( 1 ) here 


ey  = 


1 


K,r 


1/2 


yJTk  E 

which  differentiated  becomes 

1 


dc  = 


•v/iUr  E 


r-l/2  dK  ~ y r-3/2  dr 


( 1) 


(32) 


where  r is  the  distance  directly  ahead  of  the  crack  measured 
from  the  crack  tip.  The  strain  differential  varies  with  the 
position  of  sampling  r.  In  this  model  we  calculate  the 
value  of  de  at  the  fixed  point  r = d^ , and  assume  that 
strain  is  constant  from  this  point  back  to  the  actual  crack 
tip.  The  effect  on  the  local  strain  of  moving  out  from  d^ 
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a distance  dr  is  taken  to  be  equivalent  to  holding  a poin  + 
fast  in  specimen  coordinates  and  allowing  the  crack  to 
advance  a distance  -da  up  to  it,  whence  dr  = -da.  With 
this,  Eq . (32)  becomes 


di 


d 


XdK' 


( 33) 


The  second  K in  Equation  (33)  carries  a suffix  G to  denote 
its  function  as  a gradient  factor.  The  other  K term,  dKp , 
is  associated  with  rising  load  application.  It  is  needed 
later  in  connection  with  a strain-rate  field  mapping. 

However  for  corrosion  fatigue  growth,  it  may  be  neglected, 
as  the  growth  event  is  executed  during  the  load  dwell,  while 
K;  is  constant.  The  variable  Kq  is  proportional  to  the 
gradient  -uL  of  the  strain  field  at  the  point  dip,  of  which 
a s pec ial  p r ovi s i on  pertains  under  conditions  of  cyclic 
loading,  as  discussed  later.  We  retain  these  distinctions 
in  revertin?  Eq.  (33)  to  K-proportional  strains  via  Eq.  (31). 


dep  = di. 


2d, 


da 


(3M 


Eliminating  de"  between  Eq.  ( 3U  ) and  ( 31 ) , 
dcr  yields 


da  = 


8VS  dT  + 2m  dT  iji 


C,:  — - 


.-1 


Expressing  (35)  in  incremental  form  yields 

X 


Aa  = 


8V.AI  + 2m  d 


«„  (1  * 


neglecting 


( 35) 


(36) 


As  noted  earlier,  crack  propagation  per  cycle  Aa(=Aa/AN)  is 
displayed  vs  the  stress  intensity  factor  excursion,  which 
in  this  model,  from  Eq.  (l)  is 


AK  — ^/  2 7T  dp  LAc  k 


( 37) 


Here  the  0-subscript  K denotes  strain  amplitude  derived 
from  the  K amplitude.  Conversion  of  AK  or  A e (excursions) 
to  maximum  values  is  defined  in  terms  of  the  load  ratio  R 


K 


max 


AK 

1 - R 


Ac 

emax  J — 


(38) 


All  of  the  variables  in  the  right-hand  bracket  of 
Eq.  (36)(as  well  as  m)  are  measured  from  the  ordinary  or 
the  cyclic  stress- strain  curves  of  the  material,  in  a 
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manner  to  be  described  later.  The  measured  tensile  values 
are  converted  to  true  stress  o'  and  true  longitudinal  strain 
TTj  for  insertion  in  Eqs.  (36)  and  (37)-  However,  the  form 
of  the  cyclic  stress-strain  curve,  with  Bauschinger  effect, 
is  smoothed  out,  rendering  it  markedly  different,  less 
angular,  than  that  for  the  virgin  material. 


A . Dry  FCP  Growth  Extreme 


We  associate  the  equilibrium  cyclic  curve  with  fati- 
gue crack  propagation  in  a benign  or  vacuum  environment, 
applying  the  stress  relaxation  term  of  Eq.  (36).  If  we 
index  the  cyclic  flow  properties  by  numeral  2_,  as  likewise 
the  growth  associated  with  it,  the  appropriate  parts  of 
Eq . ( 36 ) are 


Ail  2 


2m  dT  tn 


( 39) 


where  £q2’  a r>  ■>  and.  0 g are  functions  of  the  true  tensile 
strain  of  the  cyclic  stress- strain  envelope.  It  is 
notable  here  that  time  enters  only  in  a relative  sense, 
a function  of  the  form  of  the  loading  pattern,  not  its 
frequency,  as  is  consistent  with  the  frequency  independence 
of  fatigue  crack  propagation  in  vacuum.  The  frequency 
independence  of  air-environment  growth  is  related  to  this 
but  has  further  implications  to  be  discussed  later. 


B . Sustained  Load  Cracking  Extreme 


At  the  other  extreme  we  associate  the  ordinary 
virgin-material,  first  (quarter)  cycle  properties,  denoted 
by  subscript  numeral  1^,  with  the  extreme  case  of  steady 
loading,  static  fatigue  under  severe  environmental  attack. 
Here,  only  the  Vg  term  is  applicable  since  with  its  linear 
dependence  on  time  it  is  generally  much  larger  than  the  m 
term,  which  decreases  with  time.  The  appropriate  part  of 
Eq . ( 36  ) is 

r / — \~i  ~ ' 


•Ail  | 


8VsAt 


pi 


(UO) 


where  egq,  Hi,  and  Oq  are  functions  of  the  true  tensile 
strain  Tq  of  the  ordinary  stress  strain  curve.  The 
coefficient  /T/ 2,  vs  unity  of  Eq.  (39)  is  to  correct  the 
instability  condition  for  plane  strain  triaxiality,  as 
proposed  by  Williams  and  Turner  [ 1 T ] • It  has  the  effect  of 
deferring  the  instability  point  to  a later  strain,  for  which 
Opl/crq  = 0.866.  The  data  correlations  indicate  that  the 
/T/ 2 coefficient  is  unnecessary,  if  not  undesirable,  for 
cyclic  growth  predictions.  Hence,  it  is  reasoned  that  the 
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reverse  straining  due  to  cyclic  loading  has  the  effect  of 
"limbering  up"  the  ligament  so  as  to  relax  the  local  tri- 
axiality.  The  /3/2  coefficient  is  retained  for  analysis  of 
stress  corrosion  cracking  (in  Eq . U 0 ) because  for  certain 
high-strength  raaraging  steels,  which  show  early  maximum 
load  points  followed  by  a protracted  flat-topped  stress- 
strain  curve,  it  defers  the  instability  point  beyond  the 
maximum  load  point  to  strain  levels  consistent  with  a 
fracture  toughness  and  process  zone  size  indicated  by 
fatigue  crack  growth  modelling. 

C . Triaxiality  Effects 


Consistent  with  the  instability  triaxiality  cor- 
rection at  the  Kjc  limit,  the  yield  condition  is  also 
assumed  to  be  affected  by  local  triaxiality.  A Tresca 
yield  condition  applied  to  the  elastic,  plane-strain  con- 
ditions which  should  border  the  plastic  zone,  equal  in 
size  to  drj,  at  this  point,  is  ; ax  = 0^.;  oz  = 2v0y.  This 
gives  a difference  of  principal  stresses  Oy  - az  = a y(l-2v) 
or  a triaxiality  factor  of  l/l-2v.  This  is  introduced  by 
augmenting  the  measured  elastic  nlus  plastic  strain  by 
( 1 _ \ 'Iys  2 v °\s 

\1  - 2v  ‘j  E 1 - 2v  E 

where  0yS  is  tensile  yield  strength  and  v is  the  elastic 
Poisson  ratio.  In  Equation  (37)  then,  Ae^-  is  augmented 


. , 2v  (Ws 

“fcK  - eKi  - ci  1 - 2v  E 


(Hi) 


and  in  Eq.  (Ho)  eG1  = eK1-  This  triaxiality  factor  for  the 
yield  condition  has  the  effect  of  elevating  the  K-propor- 
tional  strain,  via  Eq.  (37),  at  which  p_lastic  flow,  hence 
non-infinite  plastic  strain  hardening  0pi > commences.  An 
infinite  means  infinite  stability,  hence  zero  growth  by 

Eq.  (Ho).  Numerous  cases  of  severe  corrosion-induced  crack 
propagation  appear  to  be  correctly  predicted  by  use  of  this 
threshold-offset  formulation.  However  it  does  not  always 
work,  generally  underestimating  the  Kjscc  threshold  for 
softer /tougher  materials,  one  of  the  deficiencies  for  which 
this  study  suggests  remedy. 

6.  OPTIONS  IN  MIXING  OF  GROWTH  RATE  FACTORS 


An  obvious  approach  to  predicting  corrosion  fatigue 
crack  growth  is  by  mixing  additively  the  transfer  functions 
which  seem  to  work  for  each  limit:  vacuum  fatigue  at  the 

lower  end  of  the  growth  rate  scale,  Eq.  (39);  and  stress 
corrosion  cracking,  Eq.  (Ho),  at  the  upper.  But  this  does 
not  work!  In  seeking  a successful  scheme  for  combining  the 
ingredients  for  growth,  a large  number  of  different  sets  of 
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modelling  assumptions  were  tried.  Briefly,  as  it  turns  out, 
there  appear  to  be  two  distinct  rate-limiting  criteria  in 
corrosion-fatigue  crack  propagation.  The  obvious  one  is  the 
surface  ann ih i lat io n / c or r o s i on  rate  Vg  in  the  general  con- 
text of  Eq.  (Uo).  Here  the  correct  selection  of  crack  tip 
condition  is  critical.  We  have  been  able  to  do  this  by 
exercising  only  binary,  either/or  decision  options:  either 
crack  tip  conditions  pertaining  to  cyclic  attributes;  or 
else  those  pertaining  to  uncycled,  sustained-load  attributes, 
may  be  employed.  The  second  criterion,  dependent  upon  correct 
assessment  of  the  Surface  attack  criterion  , is  one  of 
constant  ratio  of  plastic  strain  rate  at  the  crack  tip  to 
the  surface  annihilation  rate.  The  first  criterion  is 
apparently  governed  by  the  corrosion  process:  it  is  frequency 

dependent  and  sets  the  level  of  the  Stage  II  growth 
corrosion  fatigue.  The  second  criterion  is  thought  to  re- 
flect a dearth  of  fresh  clean  surface  for  environmental 
attack.  Here  the  plastic  straining,  exposing  nascent  surface 
is  outpaced  by  the  corrosion  rate,  so  that  corrosion  products 
tend  to  seal  off  the  surface,  effectively  arresting  that 
straining.  To  explain  this,  the  following  sections  will 
treat  the  manner  of  estimating  the  Vg  limit,  then  the  strain 
rate  field,  then  their  combination  leading  to  the  second 
criterion . 

7.  SPECIAL  CONDITIONS  OF  FATIGUE  CRACK  PROPAGATION 

We  have  described  rather  completely  the  terms  of  the 
sustained  load  cracking  model,  summarized  by  Eq.  (Uo). 
Conditions  surrounding  the  "dry"  fatigue  growth  are  somewhat 
more  involved,  since  the  effects  of  load,  hence* strain- 
reversal,  on  the  mean  stress  location,  and  on  the  effective 
strain  gradient  require  special  definition.  Consider 

these  now  in  turn. 

A . Mean  Stress  Centroid  Relocation 

Consider  a tensile  test.  Fig.  2 in  which  the  strain 
is  reversed  and  then  restored  repeatedly  at  various  points. 
If  the  reversal  is  attempted  just  beyond  the  proportional 
limit,  point  (a),  the  plastic  yielding  will  induce  stress 
relaxation.  Since  the  true  partial  of  stress  relaxation  is 
along  the  zero  plastic  strain,  or  elastic-load  line,  the 
cyclic  stress-strain  envelope  will  tend  to  slide  back  across 
the  abscissa.  In  a low  cycle  fatigue  test,  in  which  only 
the  strain  excursion,  not  its  absolute  position,  is  con- 
trolled, the  centroid  of  area  of  the  cyclic  stress-strain 
envelope  will  slip  down  to  rest  upon  the  zero-stress  axis 
(30).  This  effect  is  attributed  to  the  inequality  of  stress 
relaxation  rates  due  to  the  imbalance  of  tensile  and  com- 
pressive stresses.  It  is  consistent  with  the  assumption. 
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■»  ( 


Sec  . . D . and  Eq . (28)  , of  a constant  relative  strain- r at e 

sensitivity,  for  which  the  stress  relaxation  A a is 

in 

proportional  to  the  level  of  flow  stress  a . 


How  does  centroid  relocation  occur  at  the  crack  tip? 

Here  the  most  that  can  be  assumed  is  that  upon  unloading 
the  crack,  collapse  of  the  elastic  stress  field  outside  th<- 
plastic  zone  will  drive  the  crack  tip  strain  back  to  zero. 

It  may  not  be  supposed  to  drive  it  into  the  compressive 
strain  quadrants.  Thus,  although  the  conditions  for  stress 
relaxation  are  present  at  point  (a),  there  is  no  available 
permissible  path  along  which  the  shift  can  occur.  To  make 
available  sufficient  c ompr  es  s i ve  - s t r es  s scope,  an  additi  :.al 
strain,  by  simple  geometric  construction  on  Fig.  2 approxi- 
mately equal  to  Oyg/2E,  must  be  provided.  This  is  done  by 
requiring  a value  of  AK,  with  its  strain  equivalent  via 
Eq.  (l)  of  Ae^,  greater  than  the  measured  longitudinal 
elastic  plus  plastic  strain  from  the  lower  toe  of  the  cyclic 
stress  curve.  The  amount  added  is  simply  the  constant  value 
o_f  ° ? S / 2E . Thus  for  a given  total  cyclic  strain  excursion 
Asp  (at  R = 0),  the  required  equivalent  strain  is  the.-, 

A£K2  = Ae 2 + ayg/2E. 


A further  consideration  is  required.  As  the  load  rati 
R increases,  the  need  for  an  incursion  path  into  the  for- 
bidden compressive  strain  quadrant  is  obviated.  The 
minimum  strain,  e_  . for  a given  strain  excursion  Aeo,  i : 

_ 2ni4  n • 

R Ey.  Thus  for  cycle-2  type  growth,  we  add  only  that  part 

of  / 2E  which  exceeds  . , whence 

Yo  2mm 


Ac 


Kl! 


(1  - R)Co  + 


»YS 

2E 


(U2) 


where  only  positive  values  of  the  square-bracketed  term  are 
employed.  The  consequence  of  this  in  terms  of  crack 
propagation  is  to  "freeze"  effects  of  load  ratio  on  growth 
rate  (threshold)  for  R values  above  1/3.  This  is  seen 
from  Eq.  (U2)  by  substituting  o^g/E  for  Ae0,  making  its 
term  zero  at  R = 1/3.  This  attribute  is  desired  for  vacuum 
fatigue  data,  for  which  R-effects  vanish  above  R - 1/3,  to 
be  demonstrated  later. 

The  next  question  is  how  to  set  the  peak-position  cyclic 
properties  for  R > 0.  Our  perception  of  this  is  depicted  by 
the  envelope  peaking  at  point  (c)  of  Fig.  2 for  an  R of  5/7. 
The  form  of  this  envelope,  following  the  Massing  hypothesis 
( 31 ) (which  we  corroborate),  is  as  though  the  full  cyclic 
envelope  were  traversed  from  its  negative  toe  to  a total 
strain  equal  to  that  of  the  strain  excursion.  Thus  the 
value  of  Oppused  in  Eq.  (39)  is  that  for 
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02  = 02(Ac2)  where  Ae2  = (1-R)e2,max 

Should  then  the  strep;-  level  associated  with  the  point 
(c)  be  adjusted  downward  for  centroid  neutralization?  It 
seems  reasonable  that  this  should  be  so.  However,  our  trails 
show  this  to  spoil  some  details  of  correlation.  Hence, 
seeking  reasons  why  it  should  not  we  note  that  the  path  of 
strain  reversal  must  follow  that  of  the  full  cyclic  curve, 
of  which  it  is,  in  effect,  a part:  the  path  d-e  of  Fig.  2. 
If  we  permit  the  c-d  envelope  to  skew  down,  execution  of 
the  d-e  path  would  require  a path  excursion  into  the 
"effectively  inaccessible"  zone.  Since  this  is  not  per- 
mitted, we  justify  retaining  the  c-d  envelope  pushed  up 
against  a "roof"  defined  by  the  full  cyclic  envelope. 


For  the  purpose  of  calculating  the  position  of  this 
"roof",  we  designate  the  compressive  stress  corresponding 
to  the  negative  toe  of  the  full  cyclic  stress-strain  curve 
as  acug , the  compressive  "ultimate  strength"  analogous  to 
aTUS  f°r  familiar  tensile  ultimate.  In  our  manner  of 

testing  for  the  cyclic  envelope,  this  point  is  at  zero 
strain,  hence  OQUS  is  a "true"  stress.  It  is  then  added 
to  true  stress  values  calculated  from  tensile  load  relative 
to  the  zero  load/stress  axis.  Calling  this  stress  a,,^,  we 
define  an  effective  or  averaged  stress  for  a given  e,  total 
cyclic  envelope 


o2  *■  <vus 

°2A  “ 9 


(1*3) 


This  value  a then  replaces  a in  the  right  bracket  of  Eq. 
(39).  2A 

B . Strain  Field  Gradient 

The  remaining  variable  in  the  cyclic  growth  equation 
(39)  is  the  gradient  strain  £q,,.  In  sustained  load  cracking, 
without  unloading,  this  is  clearly  to  be  associated  with  the 
maximum  strain  level.  However,  insisting  that  maximum 
strain  values  be  used  in  the  cyclic  growth  model,  Eq . (3?), 

has  resulted  in  prediction  of  a generally  unobserved 
negative  sensitivity  of  the  FCP  rate  to  load  ratio.  To 
effect  correlations,  we  have  had  to  assume  that  the  gradient 
strain  in  FCP  must  be  associated  with  the  strain  excursion 
in  that  cycle.  Rationale  for  this  expedient  is  to  imagine 
that  plastic  strain  reversal  in  any  amount  at  the  crack 
tip  has  the  effect  of  cancelling  out  the  tensile  strain 
gradient  (at  r = d^ ) . The  new  gradient  restored  upon  re- 
versal again  to  tension  should  reflect  only  the  positive 
strain  excursion  from  this  "zero  base",  not  that  from  trie 
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original  zero  strain,  zero  base.  Thus  the  K-equivalent 
strain  excursion  AeK2  is  used  to  replace  e„2  in  Eq.  39, 
which  for  R > 1/3  equals  Ae2.  The  use  of  XeK2  vs  Ae2  for 
R < 1/3,  is  again  an  expedient  in  view  of  the  Bucci 
corrosion  fatigue  data  at  high  R levels,  to  be  discussed 
later. 

C . Choice  of  Co  net  itions  for  Transitional  Growth 

Factor 

Crack  tip  conditions  for  cyclic  in-induced  growth 
are  likely  to  affect  those  for  subsequent  environmental 
augmentation.  The  continuation  of  the  m-induced  growth  must 
surely  be  at  the  same  local  stress  level  anc*  the 

same  local  strain  gradient  Ae^n.  However,  in  estimating  the 
centroid  neutralization  displacement,  the  first  cycle  curve. 
Fig.  2,  always  has  an  origin  of  reference  at  the  zero  strain 
ordinate.  It  cannot  move  out  from  here  as  does  the  cyclic 
envelope,  whose  shape  is  strain-excursion,  rather  than 
maximum-strain,  dependent.  Thus  the  calculation  of  Ae^., 
from  Ae^  omits  a minimum  strain  adjustment  (R/l-R)  Ae, . 
Otherwise  its  calculation  is  identical  to  that  for  Ae^ - . 
Specifically^  from  the  first  cycle  true  stress  o ^ , vs 


true  strain  £p  curve  we  define 

Aej  = (1  - FDeq 


AeKi  - ^el  2E 


This  value  of  Ae  is  converted  to  AK  by  its  insertion  in 
Eq.  (37),  to  be  aligned  (by  interpolation)  with  values 
obtained  by  inserting  Ae  from  Eq . (U2)  for  the  purpose  of 

summing  growth  rate  effects. 

Obviously  all  transition  factors  cannot  be  the  same, 
or  the  predicted  CFCP  characteristics  would  not  differ  from 
those  of  FCP.  The  difference  is  due  to  the  strain  hardening 
rate,  for_which  first  cycle  properties  pertain.  On  the 
curve  of  vs  c^,  the  0 values  are  read  as  a function  of 
e , the  maximum  strain  value.  The  true  plastic  strain 
hardening  Op,  is  then  computed  from  0^  , according  to  formulas 
given  in  section  10. 


The  environmental  attack  effective  in  one  cycle  is 
hard  to  visualize  as  a constant  corrosion  rate,  Vo. 

the  total  surface  annihilation  depth  averaged 
over  the  cycle  time  is  used  in  the  transitional  growth 


20 


formulation,  A being  substituted  for  V.At  of  Eq . 1*0. 
8.  SUMMARY  OF  GROWTH  FACTOR  FORMULAS 


Preceding  sections  cite  ingredients  of  the  three 
different  "growth  rate  factors"  now  collected.  Assume  we 
have  measured  ordinary  tensile,  then  full  cyclic-equilibr ium 
s t r es s -s t r ai n curves  and  have  converted  their  measurements 
to  "true"  values  of  o’  and  0^  as  functions  of  true 
longitudinal  strain  e.  Let  subscript  1 denote  first 
excursion  properties,  2 denote  cyclic  properties  and  2 / 1 
those  selected  as  germane  to  the  transitional  stage  between 
dry  fatigue  2 and  constant  load  cracking  1.  The  specific 
equations  are  cited  below  in  the  order  in  which  they  govern: 
dry  cyclic  fatigue,  cyclic  fatigue/corrosion  augmentation, 
and  the  static  fatigue  corrosion  induced  crack  growth.  The 
bracketed  strains  indicate  the  functional  dependency. 


For  dry  FCP,  Eq.  (39)  becomes 


Ha  2 = 2m  dT  in 


-^CK2 


Ae0  ) \ 

1 


' 2 A 


(<2) 


-1 


and  Eq . ( 37 ) 


AK  = ^/27t  dT  EAeK2 


(1*6) 


(1*7) 


in  which  , from  Eq.  1*2 

Acj^2  = (1  — R)c2  + 


«vs 

2E 


Rio 


(1*8) 


where  only  positive  values  of  the  square  bracketed  term  are 
employed . 


Fcr  the  cyclic  fatigue  corrosion  augmentation,  Eq. 
modified  as  discussed  in  Sec  7.C.  becomes 


(8ArT) 


V Ci> 
»2a(c1  ) 


(1*0) 

(1*9) 


and  Eq.  (37) 


AK  = \/2n  dT  EAeKj 


( 50) 
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where 


Ai  K , (I  K ) 1 1 


”vs 

2K 


(51) 


In  Squat. 3 ■ 1 , 

; n fa  ;1  a second 
encing  of  o 2 A ( T 


°2A  1 

excursio n 


; 

LI. 

growth , Eq . 

v" 

At 

i.  (37 


given  as  a function  of  e , when  it 
(cyclic)  property.  The  equival 
) with  the  appropriate  7^pa(e2^  is  carrle<l 

described  in 
induced  crack 


an  Interpolation  scheme  to  be 
For  static  fatigue  corrosion 
( U 0 ) becomes 

n-> 

yft 

2 


8Y 


K1 


K — TTp  Kc 


K 1 


(52) 


(53) 


where  repeating  Eq . (i+l) 


K 1 


2v  °YS 
c>  1 -2k  K 


(5M 


For  convenience  in  further  development,  that  portion  of 
each  growth  equa t ion  derived  from  the  form  of  the  stress- 


strain  curve  (except  for  m)  is 
G,  specifically: 


A . Procedure  for  Mapping 
Curves 


called  a growth  rate  factor 

(55) 

(56) 
( 57) 

Par am e trie  Growth  Factor 


The  parametric  curves  of  the  frequency-dependent 
part  of  CFCP  are  obtained  by  adding  to  the  basic  growth 
A a,.., , Eq.  (U6)  or  actually  a function  directly  proportional 
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to  it,  Eq.  ( 55  )* successively  regularly  larger  amounts  of 
Aa0  / ^ , Eq . (U9),  or  its  proportional  function,  Eq' 

(56).  The  sum  required  is: 


Aa 

(m 


•Ail  •) 


^a2/l 


2m  dT  fcn 


G 


2 


+ BArTG2/i 


(58) 


l2(m)G2  f 1 1 ( Arp ) Ci 2 / 1 


( 59) 


where  Aa  is  the  predicted  crack  growth  increment  for  one 
fatigue  cycle.  Since  G ^ is  basic,  we  factor  out  its  c o e f f i - 
c i ent 


Aa 

f2 


= G 


' 2/ 1 


(60) 


In  mapping  these  parametric  curves, 
series  for  the  coefficient  of  G 0 ^ , 

1 1 

— = 2N  X 1(T3  0 < N < 15 

1 2 


we  substitute  a geometric 
specifically 


(61) 


This  series  of  values  for  the  coefficient  of  Gp  results 
in  a family  of  curves  which  span  the  range  of  observed  CFCP 
behavior  for  the  materials  considered.  For  future 
reference  a particular  member  of  this  family  will  be  design- 
ated 


Aa 

f2(  m) 


G(N) 


(62) 


The  actual  numerical  summations  of  G(N)  must  be  made  at 
equal  values  of  A and  A c „ ^ ^ i 3 is  accomplished  by  an 

interpolation  to  be  detai lea" later . 

9.  STRAIN  RATE  FOR  CRACK  PROPAGATION 

The  concept  of  such  a rate  of  straining  inheres  this 
modelling  logic,  as  it  is  the  means  by  which  the  stabilizing 
strain  is  provided  from  crack  growth.  Recall  from  Eq.  3^  , 
crack-tip  strain  rate  combines  effects  of  crack  loading  and 
crack  growth  as  follows 
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dc  + 2d^  da 


noted  earlier,  the  gradient  strain  £q  for  CFCP  is  taken 
to  be  Ae jf . Substituting  also  the  crack  propagation  incre- 
ment in  the  context  of  Eq.  (62),  Eq . (3*0  becomes,  in 

incremental  form. 


Atp  - AcL 


Mm)  G(N) 


We  have  observed  that  the  map  of  plastic  strain  rate  rather 
than  total  elastic  plus  plastic  strain  rate  provides  the 
more  "data-oons i - lent"  delineation  of  the  frequency 
independent  contour  lines.  Accordingly  we  seek  to  remove 
the  elastic  strain  component.  To  do  this  we  deduct 
/ E from  the  cyclic-envelope  elastic  strain 


= Ac  2 


2 AP 


- f At  K m tn  I l + — G(N) 


when  ■'  ■ Inserted  ’rom  (58)  and  (5°)  and 

cf2A  is  a function  of  the  strain  excursion  Aep  n°t  the 
maximum  strain,  as  used  in  connection  with  Eq.  ( U 9 ) . In 
order  to  compare  with  a map  of  . ' , Eq.  (64)  can  be 

rearranged 


min  1 + — 


Acp  “ A<?2  v 


G(N) 


bbr<  ri ate  the  Left  and  terr  as  if.  Note  further  that  AT 
can  be  expressed  as  a strain  rate  by  normalizing  over  the 
cycle  duration  At 


"^fp 


Ac., 

P T 

At 


where  t is  the  period  of  a cycle  of  frequency  f.  With 
this,  end  using  Eq.  (U2).  Eq.  65  becomes 


m In  I + 


G(N)  h V = 


Vf  _ At2  + 2ii2A/K 


Ac2  + oYs^2K  - Ac2  {>0  only} 


This  is  shown  as  step  (19)  in  Table  1. 
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A . Procedure  tV-r  Mum  inr  i a m r i ■ • i'-i  i !.  iu'  ■■ Curves 

To  iijap  out  Eq.  (67)  we  establish  a parametric  set  of 
values  for  ep/f  using  the  same  type  of  geometric  series  as 
for  the  G t N ) map.  The  values  used  are 


_ 2m 
Vf  ~ To 


-10  < M < 4 


(68) 


for  a map  area  sufficient  to  cover  most  ranges  of  CFCP 
growt  h . 


The  maps  of  growth  rate  factor,  G(N)  vs  At,.,  and 
those  of  constant  crack  strain  rate,  t / f vs  Ac,.,  can  be 
superposed  by  aligning  coordinates.  ?f>e  abscissa  scales, 
both  in  terms  of  , mav  be  superposed.  For  the  ordinate 

scales,  we  align  Y = 1 &og  Y = 6)  on  the  strain  rate  map  in 
accord  with  Eq.  (65) 


G(N)  = 


1 


m l!n 


(69) 


Alternatively  the  strain  rate  maps  may  be  plotted  directly 
in  terms  of  G(N)  by  inserting  known  values  of  m and  t. 


B . M at  ch i ng  of  Comb i ned  Cu ryes  to  CFCP  Pat  a 

Superposition  of  the  combined  map  on  CFCP  data  plot 
produces  a value  of  dT  by  the  location  thereon.  The  pro- 
cedure used  is  to  solve  the  growth  equations  (62)  and  (hj) 
for  a line  on  the  G vs  Ae^;  map  of  constant  dp.  Since 
G var  i es(  inversely  Vith  dip  while  Ae^  as  its  root,  this 
line  has  a slope  +2  on  equally  scaled  logarithmic  paper. 

Our  preference  is  to  use  graph  paper  with  an  abscissa  scale 
twice  that  of  the  ordinate,  on  which  this  slope  appears  at  a 

1+5°  angle.  If  we  select  a reference  point  Aal  _ pg  / 

AT  ref 

cycle,  AK|ref  = 10  MPa  /m  on  the  data  plot,  a correspond- 
ing point  on  the  G vs  Ae^  match  line  is,  from  substitution 
in  Eqs.  (62)  and  (U7)  and  rearranging  with  dp  set  at  100  pm 

GU  = j— : (TO) 

20m  Sn  (l  ♦ 

AeK'-f=7^17  (inMPa)  (Tl) 
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4at  tl  a best  fit  with  the 

match  line  of  slope  2 thr  igh  Ini  ( A f v.  , G)  . r„  .• 

pas:  Lng  thr  igh  the  (&K,  Aa/AN)|r(  : Lnt.  'Using  this 

specific  value  of  A e ^ i r e f (any  other  could  be  chosen)  is 
convenient  in  that  in  match  location  it  corresponds  on  the 
AK  scale  to  /dtp  (in  pm).  The  match  must  also  honor  the 
location  of  K=  Kjq  which  should  correspond  to  the  AK/(l-R) 
at  which  G- 

10.  PROCEDURES  FOR  OBTAINING  PLASTIC  FLOW  DATi 

The  manner  of  processing  data  to  obtain  the  parametric 
growth  rate  factor  maps  will  depend  on  the  form  of  the  data. 
Hence  we  interject  at  this  point  a brief  description  of  how 
ours  was  obtained.  Tensile  and  full-cycle  equilibrium 
tension-compression  stress-strain  curves  of  the  material  are 
needed.  In  the  data  shown  here,  the  cyclic  curve  is 
obtained  after  an  initial  tensile  excursion  to  a strain 
slightly  exceeding  that  for  maximum  load.  This  allows 
measurement  of  0^  to  the  slightly  negative  values  required 
to  make  0pl/^q,  reach  the  value  /T  /2  in  Eq.  (57).  Sub- 
sequent to  this,  straining  is  reversed  and  initial  length 
is  restored.  At  this  point  the  cyclic  measurement 
commences.  Again  the  specimen  is  extended  to  slightly 
beyond  maximum  load.  In  the  ferrous  alloy  used  here, 
this  curve  is  equivalent  to  the  equilibrium  cyclic  curve. 
However,  in  the  titanium  alloys,  further  cycling  will 
typically  introduce  the  two-stage  hardening  effect.  We  had 
sought  to  avoid  this  effect  in  early  modelling  attempts  as 
there  appears  to  be  no  undulation  in  FCP  which  corresponds 
to  that  reflecting  two-stage  hardening  in  the  G (N) 
profiles,  the  parametric  lines  of  constant  "N"  in  the 
growth  rate  factor  maps.  It  is  now  apparent  that  the 
constant  strain-rate/corrosion-rate  profiles,  lines  of 
constant  N-M,  will  fit  the  data  evenwith  inclusion  of  two 
stage  hardening.  For  normal  growth  increments,  much  smaller 
than  dtp , each  ligament  should  experience  many  strain  cycles. 
Consequently,  the  condition  of  two  stage  hardening  should 
occur  at  the  crack  tip.  Accordingly,  we  now  recommend 
continued  cycling  in  these  measurements  until  a true 
equilibrium  condition  is  established.  The  strain  excursion 
applied  should  not  exceed  that  for  the  equilibrium-condition 
maximum  load,  to  assure  maximum  endurance  of  the  specimen. 

It  is  important  that  the  flow  measurements  be  made  on 
specimens  which  have  not  begun  to  fracture,  which  is  favored 
by  a minimal  strain  excursion.  The  tensile  specimen  axis  is 
oriented  normal  to  the  crack  surface,  to  simulate  dip-ligament 
orientation.  Miniature  tensile  specimens  are  employed  with 
test  section  dQ  = b . 32  mm  diameter  and  length  about 
13  mm.  Deformation  is  sensed  with  a diametral  bi-lobed  clip 
type  strain  gage,  although  we  now  believe  that  longitudinal 
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strain  measurement  would  be  more  desirable.  This  would 
remove  the  necessity  to  convert  to  longitudinal  strain  for 
use  in  the  model  calculations  providing  better  definition 
in  the  growth  threshold  region.  Usually  several  tests  of 
the  same  material  are  run  to  assure  reproducibility  of  the 
data.  The  cyclic  envelope  is  held  in  the  positive  longitu- 
dinal (=  negative  diametral)  strain  quadrants,  as  noted 
earlier,  to  better  simulate  the  crack  tip  condition  of  no 
potential  for  compressive  strain.  The  two  curves  scaled 
in  terms  of  engineering  stress  a = P/A  and  diametral  strain 
= Ad/d0  are  marked  at  closely  spaced  strain  stations, 
particularly  in  the  low  strain  and/or  high  curvature  region::  . 
At  each  station,  the  tensile  stress  Oip  and  slope  do^/de^  are 
measured  and  recorded.  Specific  point  measurements  on  each 
curve,  distinguished  as  before  by  subscript  1 for  ordinary 
and  2 for  cyclic,  are  required: 

Tensile  yield  strength,  Oy  r.  (at  0.2  $ offset); 

Elastic  diametral  moduli,  G^,-,  -y , G^q  q ; and 

Compressive  Ultimate  Strength,  Ogyg 


Additionally  a value  of  the  longitudinal  elastic  modulus  E 
is  needed;  one  may  normally  employ  handbook  values.  This 
permits  the  Poisson  ratio  to  be  calculated. 


i = 1,2 


n>.  i 


(72) 


Values  of  v specific 
those  for  the  cyclic 


°d0 , 1 * 


to  each 
toe  are 


curve  must  be  retained  since 
sometimes  higher,  i.e.,  G^q  0< 


Data  for  determining  the  stress  relaxation  exponent  m 
is  obtained  by  abruptly  stopping  the  testing  machine  head, 
with  load  recorder  gain  increased  by  a factor  of  10,  and 
observing  the  consequent  load  slackening  with  electrically 
imparted  time  marks  superposed  on  the  strain  signal.  This 
is  usually  done  in  a region  of  low  strain  hardening  rate  to 
obviate  errors  from  this  source,  noted  earlier.  A cross 
plot  of  log  stress  decrement  vs  log  time,  provides  a slope 
of  value  m in  accordance  with  Eq.  (25). 


A.  Conversion  to  True  Values 


The  load  and  diametral  contraction  data  are  con- 
verted into  true  stress  and  true  longitudinal  strain  using 
the  following  formulas,  which  can  be  found  in  references  on 
low-cycle  fatigue  techniques.  For'  clarity  in  this  section 
d and  £ represent  diametral  and  longitudinal  respectively. 
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As  previously  noted  the  l-subscr ipt  designation  is  omitted 
f ir  brevity of  other  equation  formulations.  Subscripts  e and 
p represent  elastic  and  plastic.  The  diametral  strain  is 
first  converted  to  "true"  strain  because  the  formulas  which 
follow  are  not  valid  otherwise: 


I" 


(ID 

D \d 


tn 


d0  + Ad 


(73) 


c(1  = Vn  (1  + e(1) 


(7M 


in  which  - Ah-  j which  is  the  measured  engineering  strain 
do 

The  true  stress  is  comouted  from  the  uniaxial  load 


a = P/A  - 


P A0  _ 

A0  A = °T\d) 


(75) 


o = uT(l  + ed) 


-2 


(76) 


in  which  aT  is  the  measured  engineering  stress.  The  con- 
version of  diametral  strain  to  true  longitudinal  strain  is 
effected  by  considering  the  elastic  and  clastic  components 
separately 


ed  ~ ed,< 


ed, 


a_ 

On 


Cd,p 


(77) 


in  which  0g  is  the  elastic  modulus  from  the  stress /diametral 
strain  measurement  for  the  excursion  (first  or  cyclic)  being 
considered.  Therefore 


cdp  = «d 


On 


(78) 


For  plastic  strain  only,  Poisson's  ratio  v = 0.5»  and 

e<l'P  - n—  (79) 


e?,p  = 


“ 2cH  , 
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The  longitudinal  strain  is  similarly  decomposed 


H'  _ eV,o  + e v , p " [£  2cd,p 


(80) 


(81) 


Since  the  elastic  moduli  are  related  through  the  elastic 
Poisson's  ratio,  Eq.  (72),  we  can  express  Eq.  (8l)  as 


e t = (1  “ 2i>c) 


- 9., 


cd 


(82) 


Finally,  the  ^instantaneous  slope,  0g , of  the  true  stress 
true  longitudinal  strain  curve  is  computed  using  the 
following  relationships: 


«d 


dd  _ ^fcd 
ded  ded  ded 


LU  +ed)2 


2«t 

(l+ed)3- 


(>  +cd> 


(83) 


( 8U  ) 


0a 

<1  +fd) 


2oT  _ ^d 
(l+ed)2  U+fd) 


(85) 


in  which  0^  is  the  measured  slope  of  the  engineering  stress/ 
diametral  strain  curve.  From  this  quantity,  0g  is  com- 
puted : 


r _ do  _ dd  ded  _ J_  : 
c deg  ded  deg  2 d 


(86) 


1 + 


ed 


- 2d 


29 


11.  COMPUTER  PROGRAM  FOR  GROWTH  AND  STRAIN  RATE  FACTORS 


Table  I provides  a set  of  computational  formulas  about 
which  the  computer  program*was  fashioned.  : ach  rtep  is 
numbered  so  that  the  following  comments  can  be  ea  1 ly  re- 
ferenced, although  most  steps  should  be  cel  f • vi  :•  . 

Generally,  calculations  pertaining  to  fir.t  :y  . r 

ordinary  tensile  properties  are  on  the  left  * • * at . ; 

cyclic  properties  on  the  right. 

Step  1)  converts  nominal  stress  to  true  i 
extreme  right  provides  the  centroid  normal  ir.ed  y . 1 - . tress 

°2A ' 

Step  2)  converts  to  true  strain  hardening,  s‘ill 
diametral . 

Step  3)  converts  to  true  longitudinal  elastic  plus 
plastic  strain  hardening. 

Carrying  on  through  Step  6\  two  lists  of  true  stress 
(ap,  o^,  ) , true  strain  hardening  values  (0-^,  ) versus  true, 

R-adjusted  strain  range  (Ae-,  , Aig)  are  established.  However, 
their  strain-wise  positions  do  not  necessarily  coincide,  and 
in  order  to  proceed  with  the  addition  in  G(N)  , we  must 
achieve  this  coincidence.  It  would  be  possible  to  inter- 

on  one  set  with  respect  to  the  other,  but  for  greater 
uniformity  and  accuracy  in  a subsequent  interpolation,  we 
interpolate  both  Ae-^  and  Asp  bases  with  respect  to  a third 
station  set  — a geometric  series  with  intervals  of  log- 
arithmic spacing  of  5 1 ■>  and  a strain  range  which  covers  that 
of  these  materials.  That  operation  is  given  in  Step  7),  and 
the  interpolation  is  carried  out  in  Step  8).  It  is  helpful 
to  note  that  interpolating  for  a or  0 by  addressing  the  A e 
list  results  in  values  which  correspond  to  an  e x value 
equivalent  to  Ae/(l-R).  Conversely,  if  one  interpolates 

on  the  e „ list,  values  of  a and  0 would  correspond  to 
max  ’ - 

R-adjusted  strain  range,  Ae  values. 

Step  9)  adds  the  centroid  neutralization  strain  to  the 
measured  longitudinal  (true)  strain.  However  since  this 
correction  differs  between  the  two  sets  for  R < 1/3,  it  is 
necessary  to  interpolate  again  to  realign  the  stations  at 
which  the  summations  are  to  occur.  Step  11 ) specifies  this. 
Step  12),  removal  of  the  elastic  part  of  the  strain  harden- 
ing creating  a singular  point,  is  done  here,  after  the 
interpolations,  since  linear  interpolations  suffer  in 
accuracy  when  approaching  a singularity. 

Step  13)  is  the  answer  line,  the  two  growth  factors 
Gp/i  and  Gg  needed  for  the  mapping,  which  is  specified  in 
step  l1* ) . 

* Available  from  authors  on  revest. 
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(20)  Hot  Y vs  Si  k.,.  from  eq  (9)  above  for*  l 2M  10 


■ 


Ste;  s 15)  - l8)  pertain  to  calculation  of  Gp  and 
triaxially  augmented  e^l-  Since  it  is  not  to  be  add  d to 
that  from  another  data  set,  we  simply  compute  and  plot  these 
at  the  measured-data  strain  stations.  The  strain  e at 
which  G p goes  to  infinity  should  correspond  to  K j c in  the 
data  match . 


Steps  19)  and  20)  detail  computation  of  the  strain  rate 
maps.  For  this  paper,  strain  rate  maps  were  plotted 
manually,  although  as  time  permits  we  will  computerize  these 
as  well.  Our  maps  were  computed  only  for  the  cases  R = 0 
and  R > 0.33  to  obviate  the  task  of  producing  a different 
map  for  each  R level  when  R < 0.33.  Once  automated,  strain 
rate  maps  should  be  paired  to  the  0 (N)  maps  in  a single 

computer  run  for  the  exact  value  of  R. 


12.  "DE-LUDERIZATION"  PROCEDURE 


In  this  model,  any  tensile  strain  hardening  rate  less 
than  zero,  or  true  rate  less  than  o’,  signals  an  infinite 
growth  rate  factor.  This  occurs  ostensibly  around  the  yield 
point  of  lower  strength  steels,  where  it  is  associated  with 
the  upper  yield  point  and  "delayed  yield"  phenomenon.  In 
Bucci's  s tr e s s - cor r os  ion  cracking  data  of  a highly  tempered 
U 3 U 0 steel  [32],  the  trend  in  growth  rate,  an  early  peak 
corresponding  to  the  lower  yield  or  Luder’s  strain,  was 
observed.  However,  we  have  not  seen  this  effect  in  cor- 
rosion fatigue  data  on  it  and  other  softer  steels.  We  have 
not,  in  fact,  been  able  to  reproduce  the  Bucci  observation 
in  SCC.  However,  in  the  transitional  state  associated  with 
i s likely  that  the  upper  yield  effect  Would  be 
erased  or  "unlocked"  by  the  cyclic  exercise  of  the  material. 
The  lower  yield  plateau  is  diminished  by  the  proliferation 
of  multiple  simultaneous  bands  [33]  which  such  exercise 
would  promote.  Accordingly  we  assume  that  it  is  absent  in 
the  crack  tip  ligament.  To  remove  it  from  the  measured 
curve  we  follow  a procedure  used  by  Holloman  [ 3^  ] • A plot 
is  made  of  the  log  of  true  stress  vs  strain  values  from  the 
data  sets.  A straight  line  is  fitted  through  the  points 
beyond  the  Luder's  band  region,  with  slope  measured  to 
determine  n,  the  strain  hardening  exponent.  We  measure  the 
stress  oQ  at  ?0  = ^ strain.  _Then  other  stresses  ~o  can  be 
computed  at  arbitrary  strain  z from 


n 


and  by  differentiating  Eq . (87)  the  true  strain  hardening 


na0 

0,  = ~ 
1 e 


(88) 
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Fig.  3 — A straight  line  fit  of  a logarithmic  plot  of  true  stress  vs  strain,  of 
slope  n,  provides  an  equation  for  use  in  the  Liider  band  region  of  the  me- 
dium strength  steels  of  this  study 


Fig.  4 — The  “de-Luderization”  procedure  includes  (a)  stress-strain  curve  for 
X-65  steel  with  upper  yield  point  in  first  (quarter)  cycle;  which  (b)  con- 
verted directly  to  growth  rate  factor  maps  shows  singularity;  while  (c)  con- 
verted via  equation  from  Fig.  3 is  continuous;  whence  (d)  combined  front 
from  (b)  and  rear  from  (c)  provides  G(N)  map  used 
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Fig.  5 — Combined  growth  and  strain  rate  maps  for  X-65  steel:  (a)  Growth 
Rate  Factor,  N-map,  from  Fig.  4,  (b)  parametric  set  of  lines  of  constant 
plastic  strain  rate,  M map,  in  coordinates  of  (a),  superposed  in  (c),  whence 
segmented  curves  of  constant  strain  rate/surface  attack  rate  N-M  map  are 


Plots  of  this  sort  are  shown  in  Fig.  3 for  the  three  steeds 
of  the  present  study,  in  order  of  decreasing  strength: 

4340  steel  tempered  at  1000°F;  5Ni  steel;  and  X-6 5 line- 
pipe  carbon  steel.  The  resulting  G ( N)  parametric  curve 
set  forX-65  is  shown  in  Fig.  4,  comparing  the  "unretouched" 
set  with  that  from  Eqs . (87)  and  (88)  for  an  n = 0.25.  The 
result  of  this  is  the  "hybrid"  set,  the  front  end  of  the 
normal  plot  joined  to  the  back  end  of  the  constitutive 
equation  map.  It  is  fortunate  that  for  steels,  which  are 
prone  to  delayed  yield  effects,  the  parabolic  hardening 
rule,  Eq.  (87),  generally  fits  rather  nicely  and  indeed  the 
two  sets  of  G(N)  curves  superpose  closely  beyond  the 
Ltider ' s strain  region.  "N"  values  for  the  parametric  curves 
are  noted  thereon. 

13.  COMBINED  GROWTH/STRAIN  RATE  MAPPING 

The  "de-Liideri  zed"  X-65  provides  an  example  for  com- 
bining the  two  maps.  Fig.  5,  lower  left,  shows  a para- 
metric strain  rate  curve  set  for  this  material  . 

The  lines  of  -1  slope  toward  the  top  reflect  dominance  of 
strain  rate  due  to  propagation,  Aa  control,  turning  downward 
at  low  propagation  rates  where  the  strain  rate  due  to  crack 
loading  AK  controls.  The  crowding  ofthese  lines  to  the 
right  is  produced  by  deducting  the  elastic  strain  in 
Eq.  (64).  The  "de-LUder i zed"  G(N)  map  from  Fig.  4 is 

shown  at  upper  left,  with  the  two  superposed  at  central 
right.  It  is  possible  now  to  form  a third  set  of  para- 
metric lines  by  joining  the  upper-right  to  lower-left 
trapezoidal  apexes  with  diagonals,  as  shown  in  this  figure. 

A . Curves  of  Constant  Strain  Rate/Surface  Attack 


The  significance  of  these  lines  can  be  explained  as 
follows.  They  represent  line s of  constant  M-N  (or-(N-M)). 
Recall  from  Sec.  9.  A.  Eq.  (68)  that  the  strain  rate  para- 
metric lines  represent  constant 
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while  those  for  growth  rate.  Sec.  8.  A.,  Eq.  ( 6l ) 
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Dividing  (69)  by  (90),  and  rearranging 

r-  4 00  ..  v 
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Arx 
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ill  d p ill 


1 + — 


(91 ) 


This  means  then,  that  lines  of  constant  M-N  are  of  constant 
plastic  strain  increment  relative  to  the  surface  annihi- 
lation increment,  or  alternatively,  of  constant  plastic 
strain  rate  relative  to  surface  annihilation  rate.  These 
two  manners  of  description  are  thought  to  1 ■ relevant  to 
the  two  frequency-  independent  loci  of  CFCP  data,  respect- 
ively: of  the  air  growth  curve;  and  of  the  onset  of 
Stage  II  CF  augmentation  with  respect  to  frequency.  Dis- 
cussion on  this  point  will  follow  the  case  studies. 


ll.  PROCEDURE  FOR  COMPARING  MAPS  TO  DATA 


Several  sets  of  CFCP  data  will  now  be  compared  to  the 
model  predictions.  In  each  case  data  sheet  provides  a 
tabulation  from  ordinary  and  cyclic  s t r es s -s t rai n curves. 

The  data  has  been  processed  to  obtain  true  longitudinal 
values  to  facilitate  comparison  with  data  by  non-diametral 
strain  measurements,  which  are  the  more  customary  (and  in 
fact  to  be  preferred).  Other  necessary  and/or  useful  infor- 
mation on  the  material  is  also  tabulated.  The  figures  are 
arranged  to  show  the  combined  parametric  curves  of  the  N 
series,  M series,  and  N-M  series  for  each  R-level  on  the 
left.  Then  the  nearest  fitting  in  N,  and  N-M  curves  are 
draw  on  the  CFCP  data  to  the  right.  Functional  values 
estimated  to  best  fit  the  data  are  recorded  on  the  data 
sheet  as  part  of  the  results,  which  includes  also  the 
value  of  dip  at  the  match  position.  The  characterizing 
indexes  are  reduced:  N value  to  Arp  and  N-M  to  Ac  /Ar- 
using  Eqs . (90)  and  (9l)  respectively. 

15.  MODEL  ORGANIZATION  FOR  TITANIUM-6A1-1+V 

This  alloy  has  received  much  CFCP  evaluation  because 
of  its  extensive  utilization  in  aerospace  applications, 
with  attendant  concern  for  environmental  degradation.  Thus 
we  have  a variety  of  results  to  model. 

A.  Irving  and  Beevers  , Vacuum,  R-effects,  Fig.  6 

Irving  and  Beevers  [10]  compared  CFCP  data  for  a 
mill-annealed  product  and  two  special  heat  treatments.  The 
martensitic  condition,  quenched  and  tempered,  is  only 
slightly  harder  than  the  mi 1 1 -anneal ed  material  from  Fitz- 
gerald and  Wei  used  for  our  tensile  data.  The  correlation 
Sg  ( R ) vs  A c plots  for  their  R-level  s 0.12,  0.35,  0.6l 
is  fairly  satisfactory,  and  their  observed  R- i n s e ns i t i vi ty 
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. indi  :at<  1.  >nly  the  basic  Gg  ( R ) plot  is 
the  envir  it al  effect  is  absent . 

pa  ■ id  Wei , Air  - ffects  , Pig.  7 • 

lata  ef . w a ■ an  eni  gms  f<  r the  long  timi 

n u we  sought  to  match  it  with  constant-H  curves.  At 

const  :m’  N-M  = 9 all  principal  features  seem  coherent.  G- 
wn  f ;■  R-levi  :overing  the  range  of  the  data, 
but  to  avoid  extreme  crowding  not  every  one,  as  indicated 
on  Fig.  6.  For  the  highest  load  ratio,  R = 9.9  the  9 maps 
lie  at  the  extreme  lower  limit  of  strain  rate  map  defini- 
tion. Here  the  parametric  curve  for  an  N - - 1.5  seems  to 
fil  'iio  data.  Accurate  prediction  in  this  region  requires 
extreme  sensitivity  and  accuracy  in  the  yield  point  region 
of  the  stress-strain  curves,  which  is  lacking  in  our 

xperimenta]  t ■ ■ . . neth<  ess , hi  presenl 

:tioi  i thi  ’ •);  r<  ' i.  quite  encouraging. 

Petti  , ler  and  eppner , Salt  Water  CFCP , 

Fig.  8 . 


Again  using  the  flow  data  from  the  Fitzgerald  and 
Wei  sample,  this  data  [15]  illustrates  the  leveling  of 
"Stage  II"  gr  * 4 b with  cycle  * ime , we!  n Li  and  a 
moderately  "swept -hack"  entry  profile  of  onset  of  Stage  II. 
The  crossover  of  curves  in  transition  from  Stage  II  drops 
down  from  one  to  the  other  once  the  e„/Vg  limit  corres- 
ponding to  N-M  = 15  is  exceeded.  An  N-M°of  +9 , the  same  as 
with  the  Wei  data  on  this  alloy,  appears  to  fit  the  air 
data  and  establish  an  lower-bound  envelope  for  the  Stage  II 
growth . 


1 6.  MODEL  RGAN  IS  FOR  TI-8AL-1M0-1 V 


The  material  for  our  s t r e s s - s t r ai n curves  was  obtained 
from  the  plate  used  by  Meyn.  Bucci's  specimens  were  also 
from  this  heat,  also  provided  to  him  by  Meyn. 

A.  Meyn ' s Classic  Case  of  CFCP,  Fig.  9. 

Meyn ' s work  [36]  was  first  to  show  the  now-familiar 
trend  in  CFCP  behaviors.  The  Stage  II  threshold  appears  to 
be  definied  by  a constant  N-M  = 13.5  contour,  which  differs 
from  the  subsequent  Bucci  result. 

B.  Bucci  R-effects  in  Argon  and  Salt  Water,  Fig.  10. 

Bucci's  growth  rates  [lb]  show  an  up-trend  in  Stage 
II  growth  plateaus  with  increasing  R,  an  attribute  of  the 
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model.  The  argon  data  at  both  R = 0.05  and  0.50,  all  that 
Is  published,  is  matched  by  the  same  N-M  = 10.5,  which  is 
the  same  value  as  for  Meyn ' s air  data.  The  lateral  trans- 
lation of  these  curves,  vs  the  tight  overlap  for  the  Ti - 
6A1-UV  (Fig.  7)  is  expected  if  the  N-M  for  the  gaseous- 
environment  effect  is  high,  as  appears  to  be  the  case  here. 
The  onset  of  Bucci's  stage  II  growth  appears  unlimited  by 
an  N-M  contour,  matching  instead  the  natural  contour  of  the 
4>  ( G ( N ) ) curve  for  N - 10.5.  This  data  is  extremely  critical 
of  model  calibration,  its  digestion  providing  many  insights. 
The  irregular  excursion  for  low  M sets  of  the  t^/f  maps  is 
attributed  to  inaccuracy  in  converting  diametral"  to  longi- 
tudinal strains  in  the  region  around  the  elastic  limit. 

This  is  another  reason  for  favoring  longitudinal  strain 
measurements  for  these  predictions. 

17.  MODEL  ORGANIZATION  FOR  ALUMINUM  2 219-T851 

A.  Unangst,  Shih  and  Wei  [35],  R-effects  in  Air, 

Fig.  11. 

This  is  the  only  aluminum  alloy  included  in  the 
study.  It  is  included  to  show  the  one  cam  which  we  have 
been  unable  to  model  satisfactorily.  We  have  successfully 
modelled  a number  of  other  aluminum  alloys.  The  reason  for 
this  exception  is  unknown  to  us.  The  FCP  data  is  curious 
in  that  it  appears  to  exhibit  no  definable  stage  II  region, 
to  connect  the  FCP  data  to  the  AK  Kjc  region. 

18.  MODEL  ORGANIZATION  FOR  FERROUS  ALLOTS 

We  reanalyze  by  the  present  methodology  previously 
published  NRL  data  of  Krafft  and  Smith  [6]  on  U3^0  steel 
of  four  tempers.  To  continue  the  coverage  to  lower 
strengths,  CFCP  data  on  a 5Ni  steel  and  carbon  steel  line- 
pipe  have  been  analyzed. 

A.  L3it0/lt00°F  Temper,  CFCP  data.  Fig.  12. 

A satisfactory  match,  and,  for  the  first  time,  a 
rationale  for  the  air  data  with  an  N-M  = 13  for  a lower 
growth  rate  limit. 

B.  l431*0/600°F  Temper  CFCP,  Fig.  13. 

Here  again,  a reasonable  fit.  The  10“~  torr  data 
of  Gallagher  [37],  inserted,  is  not  quite  as  slow  as  the 
prediction  for  this  material,  although  it  approaches  it. 


C.  L 3^0/ 800°F  Temper  CFCP,  Fig.  1 14 . 


Now  with  air  effects  less  severe,  the  (N-M  = 11.5) 
air  limit  is  seen  to  fashion  the  lower  limit  curve  of 
frequency  dependent  "constant-N"  growth. 

D.  ij  3^0/l000°F  Temper  CFCP,  Fig.  15. 

With  this  high  degree  of  tempering,  the  upper  yield 
point  effect  is  restored.  In  the  analysis  it  has  been  re- 
moved by  the  " de-Liider  i zat  i on " procedure  described  earlier. 

E.  CFCP  data  of  Gallagher  and  Ryder  [38]  ana  of  Cares 

and  Crooker,  [39],  on  a Steel,  Fig.  16 . 

This  material  is  rather  stable  in  salt  water,  but 
apparently  somewhat  susceptible  to  excessive  cathodic 
polarization.  Because  of  its  low  susceptibility,  relatively 
long  cycle  durations  were  required  to  show  a measurable 
effect.  " De-Liid  er  i z at  i on"  was  required  here. 

F.  X-65  Line  Pipe  CFCP  Data  of  Vosikovsky,  Fig.  17. 

The  curves,  shown  to  illustrate  the  de-Ltlder  i z at  i on 
and  strain  rate  mapping  in  Figs.  2,  3,  and  L is  overlaid 
Vosikovsky's  data  [ U 0 ] in  salt  water,  and  in  salt  water 
with  a closely  coupled  zinc.  In  the  salt  water  the  match 
is  fairly  satisfactory.  However,  with  the  zinc  coupling, 
the  air-line  profile  prediction  lies  at  higher  growth  rates 
than  measured,  which  then  is  also  higher  than  the  Vosikov- 
sky air  data,  not  shown.  Vosikovsky  employed  a triang- 
ular, rather  than  sinusoidal,  loading  waveform.  If  one 
were  to  compensate  for  this  by  reducing  the  value  of 
?.n(l  + ttf/ti)  by  a factor  of  two,  the  match  would  be  greatly 
improved.  However  this  step  is  questionable  since  FCP  data 
on  similar  steels  using  sine  wave  loading  is  consistent  with 
the  Vosikovsky  air  data.  Also  it  is  difficult  to  under- 
stand how  the  triangular  wave  form  could  have  been  main- 
tained at  the  higher  loading  frequency  of  the  air  data. 
Regarding  waveform  effects, this  model  has  no  explanation 
for  Vosikovsky's  observation  of  reduced  CFCP  rate  with 
square  wave  loading.  Here  a greater  value  of  $.n(l  + t^/t^) 
would  be  expected,  leading  to  an  increased  CFCP  rate. 

Data  of  Dawson  and  Pelloux  [8]  on  Ti-6A1-6V-2.5N  does  indeed 
show  the  expected  upward  trend. 

It  should  be  noted  here  that  it  is  not  anticipated 
that  this  model,  as  presently  constituted,  will  be  appli- 
cable to  CFCP  in  softer  alloys.  This  limitation  is  assoc- 
iated with  greater  spacing  distances  between  strength- 
controlling microconstituents,  so  that  these  distances 
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become  comparable  to  the  size  of  the  fracture  process  z ne  , 
di|. . In  mild  steel,  this  is  associated  with  control  of  frac- 
ture by  properties  of  the  matrix,  alpha-ferrite  in  thi; 
case.  Since  the  ferrite  is  softer  and  more  ductile  than  the 
overall  material  including  its  strengthening  carbides,  the 
fracture  resistance  tends  to  be  greater  than  predicted  from 
plastic  flow  measurements  on  macro-specimens  which  reflect 
the  overall  rather  than  matrix  properties  [ h 1 ] . 

This  model  has  been  set  up  with  no  correction  for  craci 
closure  effects.  It  does  appear  that  introducing  such  an 
influence  could  provide  a closer  fit  for  the  softer  steel:: 
reported  here.  This  is  a matter  for  further  study. 

19.  SUMMARY  GRAPH  FOR  "ORGANIZED"  DATA 

It  is  proposed  that  a great  deal  of  information  about 
CFCP  behaviors  can  be  summarized  in  a simple,  yet  highly 
definitive  way  once  the  three  modelling  parameters  are  det- 
ermined. The  graph  of  ArT  vs  cycle  time.  Fig.  18*  shows  the 
degree  of  frequency-dependent  cracking.  Generally  in  mate- 
rials fatigued  at  AK  levels  below  the  Kjar,c  threshold,  the 
Arrn/x  curves  should  eventually  level  out  as  Vg  -*•  0.  The 
X-65  and  5 Ni  steel,  both  of  this  description,  appear  to  be 
so  leveling.  The  value  of  the  Arn  crossover  at  the  one- 
second  cycle  ordinate  is  a measure  of  Vg,  the  surface 
annihilation  rate.  Cases  where  Vg  has  been  determined  from 
stress-corrosion  cracking  measurements  can  be  associated 
with  the  one-second  crossover  point  as  shown  on  Fig.  18  for 
Ti-8-1-1  from  Sullivan  [ U 2 ] and  for  1*3^0  from  Ref.  [6]. 
Values  of  the  size  and  strain  rat e / cor r o s i on  rate  parameters 
are  shown  in  Table  II.  It  should  be  possible  to  reconstruct 
and  extend  any  of  the  data  sets  used  here,  and  presumably 
others,  with  these  few  characterization  parameters,  plus 
appropriate  stress-strain  data. 

20.  MODEL  COMPARISON  WITH  POOK'S  AKth(R)  COLLECTION 

There  is  considerable  practical  importance  to  the 
question  of  the  load  ratio  effect  on  the  threshold  for  FCP; 
in  fact  there  is  some  question  as  to  whether  one  exists. 

The  way  this  model-attempt  has  evolved  would  suggest  its 
existence,  inasmuch  as  a truly  elastic  limit  of  the  material 
exists.  Roughly  the  model  predicts  a minimum  threshold 
value  proportional  to  Oyg/2E  at  R = 1 , this  limit  due  to 
the  centroid  neutralization  effect.  It  reaches  a value 
3oyg/2E  at  R = 0 when  the  elastic  yield  strain  has  been 
added.  If  the  value  of  drp  is  known,  the  threshold  can 
indeed  be  directly  calculated 
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Fig.  6 — The  minor  effect  of  load  ratio  R on  FCP  rate  in  vacuum  is 
compared  to  the  m-induced  G2  growth  prediction  of  the  present 
scheme  for  the  Ti-6A1-4V  noted 
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Fig.  8 — Salt-water  CFCP  data  of  Pettit  et  al  for  a Ti-6A1-4V  is  bounded  by 
frequency  independent  Aep/ArT  lines  of  N-M  indices  shown  for  Stage  II  on- 
set as  an  upper  limit  of  growth  rate  (N-M  — 15)  and  air  environment  as  a 
lower  limit  (N-M  = 9 as  in  Fig.  7 results)  whereas  crossover  is  at  constant 
ArT  for  the  cycle  (N  set). 


Fig.  9 — Salt  water  CFCP  data  of  Meyn  for  a Ti-8Al-lMo-l  V is  bounded  in 
the  manner  of  Fig.  8.  The  lower  bound  G2-only  curve  of  the  predictive 
map  is  of  the  same  trend  as  Meyn’s  FCP  rate  in  vacuum  but  fails  to  antici- 
pate the  trend  at  low  growth  AK  levels. 
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Fig.  11  — Dry  argon  FCP  data  of  Unangst  and  Wei  is  poorly  matched  by  G2 
curves,  a failure  unexpected  in  view  of  accurate  model  predictions  of  FCP 
for  many  other  aluminum  alloys 
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Fig.  12  — Fresh  (distilled)  water  and  air  data  of  NRL  on  4340  steel  of  four 
tempers  is  shown  in  the  next  four  figures.  Here  the  hardest  condition, 
205°C  tempering,  appears  nicely  bounded  N-M  = 19  for  upper  frequency- 
independent  limit,  13.5  for  the  lower  one,  with  various  constant-N  lines 
fitting  stage  II  data  in  transition  there. 
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Fig.  13  — G-map  vs  CFCP  data  for  4340  steel  as  in  Fig.  12  but  for  tempering  at 
315°C.  FCP  data  of  Gallagher  is  much  lower  than  NRL,  data  for  air  but  not 
quite  as  low  as  present  prediction  for  m-induced  growth,  the  G2-only  line. 
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Fig.  14  — G-map  vs  CFCP  data  for  4340  steel  as  in  Fig.  12  but  for  tempering  at 
427°C.  A marked  decrease  in  the  level  of  stage  II  growth  rate  is  evident  here. 
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MATERIAL : 


FE  5MI*  ICR  a HO,  V>,  (JUENCHED  AND  TEMPERED 


REFERENCE:  A-  RYDER  AND  GALLAGHER,  REF  36,  FIG-  lb 

B-  CARES  AND  CROOKER,  REF  3%  FIG-  Lb 

NOMINAL  CHEMISTRY  UT-  9CT.  {FROM  REF  3fl> 
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Fig.  16  — Two  data  sources  on  a 5-Ni  steel  are  compared  with  a G-map  for 
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MATERIAL: 


X-bS  LINE  PIPE  STEEL 


REFERENCE:  VOSIKOVSKY,  REF  40,  FIG-  17 

NOMINAL  CHEMISTRY  UT.  PCT.  {FROM  REF  40) 
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Fig.  17  — CFCP  data  on  X-65  steel  is  roughly  matched  by  G-maps  developed  in 
Figs.  3,  4,  and  5.  However,  this  strength  level  is  thought  to  violate  the  lower  limit 
of  applicability  of  this  model  as  presently  constituted  without  reckoning  for  crack 
closure  effects,  and/or  microstructural  vs  dT  size  effects. 
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Fig.  18  — Data  organization  is  proposed  in  terms  of  this  log  plot  of  ArT  vs 
cycle  duration  t,  plus  the  Aep/ArT  values  (N-M  indexes)  defining  frequency 
independent  limits  of  Table  II.  Here  the  one  second  intercept  is  compared 
with  Vs  estimates  in  stress  corrosion  cracking  experiments.  A 1:1  slope  de- 
notes constant  Vs. 
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Such  a trend  line,  normalized  with  respect  to  AKTH  at  R = 0 
is  compared  in  Fig.  19  with  the  extensive  AKc  or  A K j j 
collection  of  Pook  [1*3].  The  Pook  data  is  normalized  with 
respect  to  the  AKip^  intercept  for  R = 0 of  a regression 
line  fitted  through  his  data  points.  This  is  done  to 
reduce  the  error  of  normalizing  to  a single  measurement 
value.  Alternatively,  one  might  normalize  to  the  inter- 
cept at  some  intermediate  R value,  such  as  R = 0.33  where 
the  centroid  neutralization  effect  vanishes  in  this  model. 
Also  at  higher  R-levels,  effects  of  crack  reclosure  are 
expected  to  cease.  This  model  does  not  include  reclosure 
effects  explicitly,  except  that  due  to  the  centroid  neutral- 
ization effect.  In  general  the  model  prediction  is  not 
inconsistant  with  Pook's  data.  From  an  engineering  view- 
point, however,  it  is  discouraging  that  none  of  Pook's 
materials,  with  the  possible  exception  of  his  60-1*0  brass, 
show  a trend  favoring  that  expected  for  no  environmental 
of,fect,  the  vacuum  line  of  Fig.  19  . Is  there  no  alloy 
ondition,  heat  treatment,  etc.  , which  obtains  an  R-  in- 
sensitivity for  the  air  environment? 

21.  MEANING  OF  STRAIN  RATE/CORROSION  RATE  CONSTANCY 

The  meaning  of  a growth  rate  threshold  controlled  by  a 
fixed  ratio  of  local  strain  rate  to  surface  corrosion  rate 
is  not  altogether  clear.  The  data  base  is  limited  and  in 
need  of  further  corroboration.  However,  assuming  that 
this  will  be  forthcoming,  it  is  interesting  to  speculate 
on  possible  mechanisms. 

There  are  two  threshold  contours  to  consider.  The 
upper  one,  the  locus  of  Stage  II  growth  augmentation, 
could  well  be  a rate  controlled  process.  Here  one  might 
regard  the  strain  rate  of  the  dT-ligament  as  a means  of 
producing  fresh,  nascent  surface,  by  the  emergence  of  dis- 
location steps  at  localized  regions  at  its  surface.  The 
rate  of  supply  of  fresh  surface  via  the  strain  rate  is 
augmented  by  the  environmental  attack,  which  induces 
growth,  hence  strain  rate.  But  this  appears  to  be  a stable 
process.  There  is  a s e 1 f- s t ab i 1 i z at i on  inherent  most 
corrosion  processes.  Corrosion  produces  inert  waste  pro- 
ducts which  tend  to  seal  off  the  surface  under  attack, 
hence  slow  the  process.  Since  the  corrosion  is  self- 
generating of  its  "fuel"  in  the  sense  that  nascent  metal 
is  fuel,  one  might  expect  that  the  s e 1 f- impe danc e of  the 
supply  of  fresh  surface  might  prove  unstable  in  the 
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negative  sense.  A bit  of  surface  fouling  slows  the  attack 
rate,  hence  the  crack  growth  rate,  hence  the  strain  rate, 
hence  the  supply  of  fresh  fuel,  hence  the  attack,  etc.  ,etc.  , 
until  the  process  is  cut  off.  In  this  way,  even  though  the 
oxide  film  is  not  impervious,  still  an  absolute  cutoff 
might  be  reached  at  a particular  level  of  strain  rate 
relative  to  corrosion  rate. 

The  lower  contour  is  the  frequency  independent  air 
line.  The  parametric  line  of  constant  Ac^/Ar^  is  frequency 
independent.  Wei  [1*U]  sees  the  initial  air-oxidation  of 
clean  steel  surfaces  to  be  extremely  rapid  for  a monolayer 
or  so;  then  stop.  Presumably,  then  this  puts  a limit  on 
what  we  call  Arrp,  which  then  limits  Ae^.  The  rationale  here 
is  not  too  satisfactory  at  present. 

22.  MICROGRAPHIC  EVIDENCE  FOR  THE  dT  PROCESS  ZONE 

There  have  been  some  published  investigations  in  which 
or  from  which  a calculated  dT  from  Eq.  2 can  be  compared 
with  quantitative  microstructural  or  fractographic  measure- 
ments. The  most  reasonable  parameters  are  the  fractographic 
dimple  size  or  spacing  and  the  inclusion  size  or  spacing. 
Such  data  of  which  we  are  aware  is  listed  in  Table  III  and 
crossplotted  in  Fig.  20.  The  trend  is  positive  although 
far  from  quantitative.  Actually  th-e  departure  of  the  trend 
from  a 1:1  match,  in  Fig.  20,  is  a relatively  small  effect. 
Thus  adding  about  1 pm  to  each  process  zone  size  calcul- 
ation, as  a fixed  correct  ion,  would  shift  the  slope  to  that 
expected.  Even  if  perfect,  the  dT  fractographic  dimple 
size  correlation  would  say  nothing  of  the  numerous  non- 
dimpled  fractographic  features  associated  with  high  cycle 
fatigue  and  s t r e s s - c or r o s i on  cracking.  It  can  be  argued 
[ 1+ 5 ] that  the  instability  event,  which  controls  the  rate  of 
the  fracturing  process,  could  be  followed  by  a separation 
process  which  tends  to  obliterate  the  evidence  for  such  a 
stage.  At  the  present  time,  the  more  definitive  way  of 
determining  the  drp  quality  index  is  by  data  vs  G-map 
matching  procedures. 

CONCLUSIONS 

The  following  are  indicated  by  this  study. 

A.  It  is  possible  to  track  CFCP  data  with  a double 
set  of  parametric  lines.  The  first  of  these  represents  the 
frequency  dependent  growth.  It  is  derived  by  summing  a 
transitional  growth  rate  factor  to  that  basic  to  vacuum 
fatigue.  The  other  is  derived  from  the  first,  overlaid 
by  a map  of  the  strain  rate  at  the  crack  tip  due  to  crack 
loading  plus  crack  propagation.  This  provides  loci  of 
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Fig.  20  — Comparison  of  calculated  dT  process  zone  size  with  measured  fracto- 
graphic/microstructural  features,  from  sources  noted  (see  Table  III)  shows  a rea- 
sonable correspondence 
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constant  ratio  of  strain  rate  to  surface  attack  rate. 

These  loci  match  CFCP  trends  in  air,  and  the  threshold  of 
Stage  II  growth  in  aqueous  environments. 

B.  A minimum-parameter  characterization  scheme  results 
from  these  considerations.  It  involves  finding  the  members 
of  the  two  parametric  families  best  fitting  a particular 
data  set.  Such  characterization  appears  to  provide  a 
useful  predictive  capability  for  estimating  effects  of 
CFCP. 
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APPENDIX 

ILLUSTRATED  SUMMARY  OF  NR L MEMO  REPORT  #3505 

1.  Corrosion  fatigue  crack  propagation,  CFCP,  is  frequently 
implicated  in  mechanical  failures  of  machinery  and  struc- 
tures subjected  to  frequent  loading  fluctuations  in  opera- 
tional environments . However,  CFCP  behaviors  vary  so 
markedly  with  material,  environment,  load  ratio  and  fre- 
quency that  engineering  estimates  of  its  course  are  impre- 
cise, and  data  required  to  refine  predictions  costly.  This 
report  describes  improvements  of  a model  relating  CFCP 

to  Engineering  mechanical  properties  - ordinary  and  cyclic 
stress-strain  curves  - by  which  many  hitherto  anomolous 
behaviors  are  reconciled.  The  new  method  could  improve 
accuracy  and  reduce  cost  of  safe-life  predictions  for 
engineering  structures. 

2.  Developed  here  is  a connection  between  the  strain  rate 
of  material  at  the  crack  tip  and  CFCP.  It  is  found  that  the 
plastic  strain  rate  in  the  loading  cycle  time  must  exceed  a 
level  c ommen s ur at e with  the  total  corrosive  activity  in  the 
cycle  before  the  environment  can  affect  crack  growth.  Thus 
the  threshold  of  lorrosion-fatigue  augmentation,  generally 

a large  jump  in  crack  growth  rate,  hence  decrease  in  "safe 
life",  is  defined  by  a locus  of  constant  ratio  of  crack-tip 
strain  rate  to  crack-tip  corrosion  rate.  Parametric  maps 
of  such  loci,  readily  prepared  from  s t r es s - s t r ai n data, 
provide  a useful  means  for  characterizing  CFCP  behavior. 

3.  Briefly  this  summary  reviews  concepts  of  the  earlier, 
basic  model,  which  provides  crack-growth  rate-factors  for 
the  limiting  cases  of  dry  vacuum  fatigue  and  of  stress 
corrosion  cracking.  A hybrid  growth-rate-factor  is  then 
developed  and  combined  with  gr owt h- f ac t or  parametric 
mapping  to  define  the  frequency  independent  growth  rate 
limit  of  CFCP.  An  example  of  model-to-data  comparison 
illustrates  the  manner  of  presentation  of  the  numerous  case 
studies  of  the  full  report. 
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U.  The  hypothesis  of  this  model  is  that  cracks  grow  to, 
and  only  to,  that  extent  required  to  maintain  a constant 
load  condition  within  a small  fixed  distance  from  a crack 
tip.  Examining  this  proposal  leads  to  two  basic  questions: 
firstly,  how  is  the  constant-load  condition  maintained  in 
ordinary  uniaxial  tensile  deformation;  and  secondly,  how 
can  crack  growth  serve  to  maintain  it  in  the  context  of  the 
corrosion  fatigue  environment?  Consider  now  these  ques- 
tions in  turn. 


Fig.  A-l  — Load  partials 

5.  In  the  ordinary  tension  test  of  a material,  free  of  an 
upper  yield  point  effect,  constant  load  occurs  in  the  max- 
imum or  ultimate  tensile  stress  region.  Here,  Fig.  A1 , the 
strain  hardening  rate,  as  a strength  partial,  has  decreased 
to  a level  too  low  to  compensate  for  load  loss  due  to 
lateral  (Poisson)  contraction,  as  an  areal  partial.  Two 
additional  partials  can  affect  crack  tip  stability:  a 

strength  partial,  stress  relaxation;  and  an  areal  partial, 
corrosion-like  surface  attack.  Three  of  the  four  partials 
are  negative,  hence  reduce  the  load  capability;  only  that 
of  strain  hardening  is  positive.  Since  strain-hardening 
requires  straining,  setting  the  sums  of  the  four  partials 
to  zero  indicates  how  much  strain  is  required  to  strike 
this  balance. 
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6.  Next  we  ask  how  such  straining  can  be  provided  to  "speci- 
men material"  within  some  fixed  distance  from  a crack  tip. 
There  are  two  ways.  Fig.  A2  depicts  that  due  to  crack  load- 
ing. Here  we  envisage  a tiny  tensile  specimen  as  connected 
to  the  strain  concentration  or  singularity  pattern  associated 
with  the  loaded  crack.  As  load  increases,  so  does  the  whole 
pattern  of  strain  concentration,  hence  so  does  the  strain. 


€ ^ 


Fig.  A-3  — Strain  from  crack  growth 


7.  Alternatively  suppose  load  on  the  cracked  body  is  main- 
tained, while  the  crack  advances  towards  a fixed  point  in  the 
material.  Then, Fig.  A3,  our  "specimen"  coupled  to  the  strain 
singularity  "ramp"  is  extended  by  the  wedging  action  of  the 
advancing  ramp.  In  balancing  the  load-partial  equation  for 
CFCP,  only  the  propagation  strain  of  Fig.  A3  is  applicable 
since  the  growth  is  executed  after  the  load  has  been  applied 
to  the  crack.  However  the  loading-strain  term  is  needed  to 
estimate  the  total  strain,  hence  strain  rate,  averaged  over 
the  total  cycle  time. 
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Fig.  A-4  — Stress/strain  — crack  growth 


8.  We  solve  for  a crack-growth  relationship  by  combining 
the  growth  vs  local  strain  relationship,  Fig.  A3,  with  the 
load  maintenance  vs  strain  relationship , Fig.  A1 , to  obtain 
a growth  relationship.  Fig.  kh  . This  is  set  up  in  dimen- 
sionless form,  in  growth  rate  factors  G which  are  propor- 
tional to  crack  growth  per  cycle,  and  a strain  AeK,  which  is 
proportional  to  crack  loading  or  the  stress  intensity  factor 
excursion  AK.  The  growth  rate  factor  curves  are  specific 
to  a given  material  through  its  particular  stress-strain 
curve(s),  hence  its  strain  hardening  behavior.  The  ordinary 
tensile  curve,  denoted  by  1_  in  Fig.  AU  , is  typically  of 
more  acute  yield  point  than  the  cyclic  curve  2_;  so  corres- 
pondingly is  its  growth  rate  factor  G. . We  regard  G0  as 
basic,  applying  to  it  a proportionality  factor  associated 
with  the  stress-relaxation  partial.  Sustained  load  or 
stress-corrosion  cracking,  on  the  other  hand,  is  associated 
with  G-]_  . For  corrosion  fatigue,  a transitional  growth  rate 
factor  Gg/p  is  defined  which  employes  the  strain  hardening 
rate  of  tne  ordinary  1_  (true)  stress  strain  curve  but  the 
stress  of  the  cyclic  2_  curve.  Since  the  cyclic  stress 
envelope  is  lower  and  smoother,  so  correspondingly  is  Gg/l 
vs  G^  only.  We  associate  G2 /p  with  the  corrosive  partial 
only,  as  an  adjunct  to  the  stress  relaxation  partial  uti- 
lyzed in  conjunction  with  Gg. 
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9.  Three  additional  features  of  the  new  modelling  should  be 
noted;  illustrated  by  Fig.  2 of  the  full  report. 

a.  Firstly,  a centroid  neutralization  rule  is  invoked, 
whereby  the  cyclic  true  stress-strain  envelope  skews  along 
its  elastic  load  line  until  its  centroid  of  area  rests  on 
the  zero  stress  axis.  For  very  low  plastic  strain,  stress 
excursions  just  exceeding  the  yield  point,  this  is  imposs- 
ible without  introducing  compressive  strain,  which  in  turn 
is  an  impossible  effect  of  an  elastic  stress  field  collaps- 
ing upon  crack  unloading.  A slight  additional  tensile 
strain,  hence  AK-offset,  one-half  the  yield  strain,  permits 
the  centroid  neutralization  under  these  conditions.  This 
strain  offset  sets  the  minimum  value  of  AK  for  the  fatigue 
propagation  threshold,  however  high  the  load  ratio  R. 

b.  The  second  feature  pertains  to  the  disposition  of 
the  cyclic  envelope  for  R > 0.  We  assume  it  to  be  of  the 
same  form  as  the  full  cyclic  curve  of  reduced  strain 
excursion,  but  subtended  from  the  profile  of  the  full  cyclic 
stress-strain  curve  exercised  to  the  maximum  strain. 

c.  The  third  feature,  closely  related  to  the  second, 
is  that  the  strain  gradient,  the  ramp  slope  for  propagation 
in  Fig.  A3,  is  assumed  proportional  to  the  strain  excursion 
As{£,  not  the  maximum  strain,  emax. 
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Fig.  A-5  — Combined  growth  sources 


10.  One  doesn't  know,  prior  to  comparing  G plots  with  crack 
growth  rate  data,  how  much  of  the  "corrosion  partial"  is 
applicable  to  a particular  mat e r i al / c orr os i ve  environment 
combination.  How  much  involves  a fitting  parameter.  To 
provide  a range  of  selections,  we  map  out  a parametric 
family  of  G curves,  adding  to  the  basic  Gg  an  arbitrary 
initial  coefficient  times  Gg/q,  Fig.  A5.  Then  a second 
plot  is  made  with  the  arbitrary  coefficient  doubled,  then 
again  times  ^ , 8,  16 , etc.  to  2^5.  We  call  the  sum  G(N)  , 
and  in  data  comparisons  ask  what  value  of  N provides  a 
member  of  the  parametric  curve  family  best  fitting  the 
growth  rate  data  of  a given  loading  frequency. 
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Fig.  A-6  — Combined  crack-tip  strain  rate 


11.  The  next  step  is  to  map  out  a parametric  family  of 
curves  of  constant  strain  rate.  This  family,  Fig.  A6  , is  of 
the  same  geometric-series  base  (two)  and  (log)  coordinates 
as  the  G(N)  vs  Ae^  map  of  Fig.  A5 • Two  limiting  condi- 
tions are  observed.  Strain  due  to  a given  amount  of  crack 
propagation  increases  as  the  K-level,  hence  gradient  steep- 
ness, increases.  Hence  lines  of  constant  strain  rate  due  to 
propagation  slope  down  -1  in  G(N)  /Aek  space.  However, 
the  loading  strain,  independent  of  crack  growth,  is  a 
function  of  AK  alone.  A base  2 family  is  a set  of  vertical 
lines  an  equal  2X  apart  on  log  paper.  We  use  only  the 
plastic  part  of  the  loading  strain  so  that  deducting  the 
elastic  strain  compresses  this  set  to  the  right.  The  com- 
bined loading  and  propagation  strain  rates  joins  these  two 
extremes.  In  the  text,  indexing  integers  of  the  strain 
rate  set  are  designated  by  M,  those  for  the  growth  rate  set 
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Fig.  A-7  — Growth  X strain  rate  mapping 

step  is  to  combine  "M"  and  "N"  maps,  by  simply 
hem  as  in  Fig.  A7 . The  loci  of  certain  inter- 
hese  now  obtain  a special  significance.  Traces 
-M  (or  M-N)  indexes  represent  lines  of  constant 
.1  plastic  strain  to  total  surface  attack  during 
ycle:  Ac  /ArT.  Relative  to  cycle  duration 

es  - plaltic  strain  vs  surface  attack  rates: 


13.  The  fit  of  the  combined  mapping  of  Fig.  AT  to  real 
CFCP  data  is  done  by  selecting  those  factors  of  proportion- 
ality providing  best  fit  of  G/Ae^  to  the  —/UK  data.  This 
fit  requires  selecting  one  parameter-  the^fiicrostructural 
size  constant  d,p . Since  dx  affects  both  ordinate  and 
abscissa  scales,  its  selection  constrains  the  fit  to  one 
degree  of  freedom:  the  reference  point  established  on  the 
growth  X strain  rate  map.  Fig.  AT,  must  lie  on  the  refer- 
ence line  of  the  crack  growth  data  plot.  Fig.  A8  . Once  dtp 
is  obtained,  we  observe  that  locus  lines  of  constant  N-M 
match  the  two  l'requency-independent  growth  rate  limits  of 
CFCP.  The  lower  one  in  AK  (upper  in  growth  rate)  appears 
to  define  the  threshold  of  Stage  II  growth,  the  onset  of 
the  frequency-dependent  plateau  region  on  the  vs  AK 
plot,  Fig.  A8.  Generally,  a higher  index  N-M,  'iaeans  a 
lower  critical  value  of  relative  strain  rate  e /Vg,  which 
corresponds  to  a lower  and  more  steeply  sloped^  ° AK  line 
for  stage  II  transition.  When  the  stage  II  transition  is 
unfavorable,  i.e.,  of  low  Cp/V^,  so  generally  is  the  upper 
AK  limit  which  corresponds  with  the  air  environment  growth 
profile . 
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Ih . To  recapitulate,  various  anamolous  behaviors  of  CFCP 
are  fitted  by  use  of  three  characterizing  parameters,  taken 
in  conjunciton  with  tensile  test  data  of  each  material. 

The  parameters  are: 

a.  The  process  zone  size  drp  , which  locates  in  one 
degree  of  freedom,  the  map  on  the  data.  It  is  the  crack  tip 
distance  within  which  the  constant  loa,  condition  is  main- 
tained. 

b.  The  surface  attack  per  cycle  ArT , associated  with 
the  N index,  for  the  growth  rate  curve  of  best  fit  to  data 
of  a given  loading  frequency.  In  susceptible  alloys,  the 
surface  attack  rate  A rip  • f = Vg  is  typically  a constant, 
independent  of  frequency;  and 

c.  The  new  parameter,  the  relative  plastic  strain 
(rate)  As  /Arp  - £p/Vg  associated  with  the  N-M  index. 
Generally^ one  value  N-M  pertains  to  stage  II  threshold,  a 
lower  one  to  the  air  environment  growth  rate. 

14.  The  efficacy  of  this  three-parameter  characterization 
of  CFCP  is  demonstrated  by  some  dozen  examples  in  the 
report.  They  include  many  of  the  known  anomalous  behaviors 
of  CFCP.  This  success  speaks  of  the  opportunity  to: 

a.  Analyse  more  cases  of  CFCP  data  with  the  improved 
model;  and  succeeding  in  this 

b.  Convert  the  method  into  simplified  engineering  rules 
and/or  computer  programs  for  safe-life  prediction  and 
damage  tolerant  design  analysis. 
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LIST  OF  ABBREVIATIONS,  ACRONYMS  AND  SYMBOLS 


K 

- 

crack  tip  stress  intensity  factor,  after  Irwi: 

Kl 

- 

mode  I,  plane  strain  stress  intensity  factor 

AK 

_ 

stress  intensity  factor  excursion 

KIc 

Kj  for  onset  of  fast  fracture  instability  or 
plane  strain  fracture  toughness 

^Iscc 

- 

Kj  for  threshold  of  SCC 

FCP 

— 

Fatigue  Crack  Propagation 

CFCP 

- 

Corrosion  Fatigue  Crack  Propagation 

see 

— 

Stress  Corrosion  Cracking,  sometimes  called 
static-fatigue  corrosion-induced  crack  growth 

P 

- 

load  on  a specimen 

R 

- 

load  ratio,  I'm|n/Pmax,  in  fatigue 

a 

- 

crack  length 

ajj  i Aa/  AN 

increment  in  crack  length  in  one  cycle  of 
fatigue  loading 

f 

- 

frequency  of  fatigue  loading 

T- At 

- 

period  of  fatigue  loading  1/f 

tL’  tH 

— 

time  intervals  of  loading  and  load-holding  of 
a fatigue  waveform 

V 

- 

Vacuum  Environment 

A 

- 

Air 

IA 

- 

Inert  Argon 

FW 

- 

Fresh  Water 

SW 

- 

Sale  Water  (3-5%  NaCl) 

SEA 

- 

Key  West  Ocean  Water 

a , a , a 

_ 

stress,  true  stress,  stress  rate  respectively 

0 YS 

- 

yield  strength  at  0.2 l offset 

°TUS 

- 

Ultimate  Tensile  Strength 

°CUS 

- 

Ultimate  Compressive  Strength 

0 

- 

slope  of  stress-strain  curve  or  tangent 
modulus 

E 

- 

Young's  elastic  modulus 

V 

- 

Poisson  ratio 

m 

- 

stress -relaxation  exponent 

e , T,  e 

- 

strain,  true  strain,  strain  rate  respectively 

ed 

- 

diametral  strain 

eZ 

- 

longitudinal  strain 

eP 

- 

plastic  (only)  strain 

£1 ,2 

- 

strain  pertaining  to  first  excursion  (l),  or 
full-cyclic  equilibrium  (2)  strain  cycle 

S 

- 

testing  machine  stiffness 

C rp 

— 

total  compliance  of  a specimen-machine  assem- 
blage 
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£y  - stress,  strain  in  direction  normal  to  plane  of 

mode  I crack 

- distance  ahead  of  crack  or  singularity  point 
of  crack-tip  stress  field 

- a specific  value  of  r , found  as  a size  para- 
meter in  the  fracture-flow  modelling 

- radius  of  a ligament  of  diameter  d-p 

- radial  diminution  due  to  environmental  surface 
attack  or  annihilation  during  the  load  cycle 

- rate  of  surface  attack,  Ar-p/x 

- cross-sectional  area  of  dp_ ligament 

- change  in  ligament  cross-section  due  to  Poisson 
contraction 

- change  in  ligament  cross-section  due  to  surface 
annihi lat ion 

- change  in  stress  on  ligament  due  to  strain- 
hardening of  ligament  material 

- change  in  stress  on  ligament  due  to  stress 
relaxation  of  the  ligament  material 

value  of  K associated  with  the  strain  gradient 
near  the  crack  tip 

value  of  K associated  with  load  application  in 
the  fatigue  cycle 

gradient  strain,  proportional  to  at  r = dT  , 
also  proportional  to  Kq  ^r 

loading  strain,  at  r = d>p , due  to  application 
of  load  in  fatigue  cycle,  also  proportional 
to  Kl 

strain  amplitude  associated  with  maximum 
stress  intensity  factor  or  its  excursion 

plastic  strain  accruing  in  one  loading  cycle 
due  to  hoth  crack  loading  and  propagation 

- growth  rate  factor  associated  with  cyclic  (2) 
stress  strain  curve 

- growth  rate  factor  associated  with  first 
excursion  (l)  stress  strain  curve 

p - growth  rate  factor  proportional  to  the  corro- 

sion-augmentation of  cyclic  fatigue  crack 
growth 

m)  - coefficient  of  G2 

Arrr,)  - coefficient  of  G,  nr  r,„ /, 


or  Aek 


f 


G(N)  _ growth  rate  factor  found  by  adding  varying 

degrees  of  &2/1  ^2 

Y - strain  rate  factor  proportional  to  . G(N)  used 

to  display  curves  of  constant  plastic  stia^n 
rate  of  the  dT  ligament,  ep/f 

jj  _ geometric  series  exponent  for  parametric  cur'."  s 

of  constant  Ar^  in  G(N)  vs  Ae^  space 

M - geometric  series  exponent  for  parametric  curves 

of  constant  £ /f  in  Y vs  A eK  space 

CONVERSION  FACTORS 


1 MNm" 
1 MNm" 


-2  =1  MPa  = 0.1^5  ksi 

1/2 

-3/2  = i MPa  • m = 0.91  ksi 


1 um 


39  . U 


